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' Abstract 

O ■ We calculate the CP-averaged branching ratios and CP-violating asymmetries of a number of two-body 

d : 

charmless hadronic decays B® — > PP, PV, VV in the perturbative QCD (pQCD) approach to leading 

b> ■ 

.£h ! order in a s (here P and V denote light pseudo-scalar and vector mesons, respectively). The mixing- 
X 1 

$_i ' induced CP violation parameters are also calculated for these decays. We also predict the polarization 



fractions of B s — ► VV decays and find that the transverse polarizations are enhanced in some penguin 
dominated decays such as B® — > K*K*, K*p. Some of the predictions worked out here can already be 
confronted with the recently available data from the CDF collaboration on the branching ratios for the 
decays B® — > K + ir~ , B® — ► K + K~ and the CP-asymmetry in the decay B® — ► K + tt~, and are found 
to be in agreement within the current errors. A large number of predictions for the branching ratios, 
CP-asymmetries and vector-meson polarizations in B® decays, presented in this paper and compared 
with the already existing results in other theoretical frameworks, will be put to stringent experimental 
tests in forthcoming experiments at Fermilab, LHC and Super B-factories. 
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I. INTRODUCTION 



There has been remarkable progress in the study of exclusive charmless B® — > hih 2 and 
— ► hifi2 decays, where h\, hi are light pseudo-scalar and/or vector mesons. Historically, these 
decays were calculated in the so-called naive factorization approach jlj] , which was improved by in- 
cluding some perturbative QCD contributions . Currently, there are three popular theoretical 
approaches to study the dynamics of these decays, which go under the name QCD factorization 
(QCDF) 4], perturbative QCD (pQCD) fl, and soft-collinear effective theory (SCET) Q. All 
three are based on power expansion in l/m^, where m& is the 6-quark mass. Factorization of the 
hadronic matrix elements (hih2\Oi\B) , where 0% is typically a four-quark or a magnetic moment 
type operator, is shown to exist in the leading power in l/m& in a class of decays. In addition, 
these approaches take into account some contributions in the decays B — > hiho not included in 

n u n 

the earlier attempts |2|, |3| , in particular the so-called hard spectator graphs. 

Despite being embedded in the A/m^ approach, justified by both the large mass, m& = 0(5 
GeV), and a large energy release in the decay, with E^. = m#/2, these methods differ significantly 
from each other in a number of important aspects. For example, these differences pertain to 
whether one takes into account the collinear degrees of freedom only as in QCDF and SCET, or 
includes also the transverse momenta implemented using the Sudakov formalism, as followed in 
the pQCD method. Also, in pQCD, the power counting is different from the one in QCDF, which 
makes some amplitudes differ significantly in the two approaches. The other differing feature of 
pQCD and QCDF is the scale at which strong interaction effects, including the Wilson coefficients, 
are calculated. In pQCD, this scale is low, typically of order 1-2 GeV. In QCDF, typical scales 
for the Wilson coefficients are taken as O(mb), following arguments based on factorization. There 
also exist detailed differences between QCDF and SCET, despite the fact that dedicated studies 
in the context of SCET have allowed to gain a better understanding of the QCDF framework. 
These differences, though not inherent, lie in how practically the calculations are done in the two 
approaches and involve issues such as the treatment of the so-called charming penguin contribu- 
tions [3| to the decays B — > h\h 2 . These are argued to be power-suppressed in QCDF, and left 
as phenomenological parameters to be determined by data in SCET. Likewise, the treatment of 
the hard spectator contribution in these two approaches is also different. We recall that a generic 
factorization formula 8| 

{hMOi\B) = $ h2 (u) * (T I (u)F Bhl (0) + C h (t, u) * E Bhl (r, 0)) (1) 



involves the QCD form factor F Bhl (0) and an unknown, non-local form factor E Bhl (r, 0). In 
QCDF, this non-local form factor factorizes into light-cone distribution amplitudes and a jet 
function J(t,u,v), when the hard-collinear scale y/m^K is integrated out. This interpretation 
of the hard spectator contribution was not at hand in the BBNS papers 4J], but was gained 
subsequently in the SCET-analysis of the form factors jj. Amusingly, this SCET-result is not 
used in the SCET-based phenomenology in B — > hih 2 decays, for example in the works of Bauer et 
al. 1Q[ , where the use of perturbation theory at the scale ^m b A is avoided. Detailed comparisons 
of their predictions with the data for the decays of the B® (and its charge conjugate) and B ± - 
mesons have been made in the literature. We also refer to a recent critique 8| of the underlying 
theoretical assumptions in the three methods. Data from the B-factory experiments, BABAR and 
BELLE, as well as the CDF collaboration at the Tevatron, do provide some discrimination among 
them. With the advent of the LHC physics program, and the steadily improving experimental 
precision at the existing facilities, it should be possible to disentangle the underlying dynamics in 
hadronic 5-decays. 

The experimental program to study non-leptonic decays B® — > h\h 2 has also started ll|, with 
first measurements for the branching ratios B® — > K + tt~ and B® — > K + K~ made available recently 
by the CDF collaboration 12J, ll3[ at the proton-antiproton collider Tevatron. Remarkably, the 
first direct CP asymmetry involving the decay B® — > K + ir~ and its CP conjugate mode has also 
been reported by CDF [13|, which is found to be large, with A CP (B® — > K + n~) = (39 ± 15 ± 8)%. 



This is in agreement with the predictions of the pQCD approach, as we also quantify in this 
paper. With the ongoing 5-Physics program at the Tevatron, but, in particular, with the onset 
of the LHC experiments, as well as the Super B-factories being contemplated for the future, 
we expect a wealth of data involving the decays of the hitherto less studied B® meson. The 
charmless B® — > h\h 2 decays are important for the CP asymmetry studies and the determination 
of the inner angles of the unitarity triangle, in particular 7 (or </> 3 ), which has not yet been 
precisely measured. In addition, a number of charmless decays B® — > hih 2 can be related to the 
Bj — > hih 2 decays using SU(3) (or U-spin) symmetry, and hence data on these decays can be 
combined to test the underlying standard model and search for physics beyond the SM under 
less (dynamical) model-dependent conditions. Anticipating the experimental developments, many 
studies have bee,, devoted to the interesting charmless 5? M> decays. Ar^ng others, they 



include detailed estimates undertaken in the naive 
generalized factorization approach [l^ . QCDF 
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actorization framework [14j | , th e so-called 
pQCD [19J and SCET [20]. There are 
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also many studies 21] undertaken, parameterizing the various parts of the decay amplitudes using 
distinct topologies and the flavor symmetries to relate the — > h-\ hi and B\ — > h\hi decays. 
Possible New Physics effects in these decays have also been explored 22j. 

In the applications of the pQCD approach to B® — > h\h2 decays, the currently available works 
concentrate on specific decays. However, a comprehensive study of the decays B® — > h\h2, which 
have been undertaken in QCDF and SCET, to the best of our knowledge, is still lacking in pQCD. 
Our aim is to fill in this gap and provide a ready reference to the existing and forthcoming 
experiments to compare their data with the predictions in the pQCD approach. In doing this, we 
have included the current information on the CKM matrix elements, updated some input hadronic 
parameters and have calculated the decay form factors in the pQCD approach. Since these form 
factors have to be provided from outside in QCDF and SCET (such as by resorting to QCD 
sum rules), there is already a potential source of disagreement among these approaches on this 
count. However, we remark that the estimates presented here for the cases B s — > PP, PV, VV 
are rather similar to the corresponding ones in the existing literature on the light cone QCD sum 
rules. Thus, theoretical predictions presented in this work reflect the detailed assumptions about 
the dynamics endemic to pQCD, such as the effective scales which are generated by the strong 
interaction aspects of the weak non-leptonic two-body decays, setting the relative strengths of 
the various competing amplitudes in magnitudes and phases. As we work in the leading order 
(LO) in a s , there is considerable uncertainty related to the scale-dependence, which we quantify 
in the estimates of the branching ratios, CP-asymmetries and polarization fractions for the decays 
considered in this paper. Likewise, parametric uncertainties in the numerical estimates of these 
quantities resulting from other input parameters are worked out. Together, they quantify the 
theoretical imprecision in B® — > h\hi decays at the current stage in the pQCD approach. We have 
made detailed comparison of our predictions with the existing literature on the decays B® — ► h\hi 



and in some benchmark decay widths and rate asymmetries in the corresponding B®, — > h±h 
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decays. Whenever available, we have also compared our predictions with the data and found that 
they are generally compatible with each other. 

We also present numerical results for some selected ratios of the branching ratios involving 
the decays B° -> K + K~, W s -> K+ir', ~B% -> tt+tt" and ~B% -> K~tt + , which are related by 
SU(3) and U-spin symmetries. For comparison with other approaches, we also give in Table IV 
the contributions of the various topologies for these decays. The ratios worked out numerically 
are: R\ = BR( - B ^ K K \ R 2 = BR ( B S ~^ K K ) anc [ ^ wo m0 re, called and A, defined in Eq. (72) 
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and (73), respectively, which involve the decays B® — > K~ir + and B® — > K + n~ and their charge 
conjugates. All these ratios have been measured experimentally and the pQCD-based estimates 
presented here are in agreement with the data, except possibly the ratio R\ which turns out too 
small. Whether this reflects an intrinsic limitation of the pQCD approach or the inadequacy of 
the LO framework remains to be seen, as complete next-to-leading order (NLO) calculations of all 
the relevant pieces of the B — > hxh-i decay matrix elements are still not in place. However, there 
are sound theoretical arguments why the ratios R2, R3 and A are protected against higher order 
QCD corrections, such as the charge conjugation invariance of the strong interactions (for R3 and 
A), and the dominance of the decay amplitudes in the numerator and denominator in the ratio 
R2 by a single decay topology. The agreement between the pQCD approach and data in these 
quantities is, therefore, both non-trivial and encouraging. 

This paper is organized as follows: In section II, we briefly review the pQCD approach and give 
the essential input quantities that enter the pQCD approach, including the operator basis used 
subsequently and the numerical values of the Wilson coefficients together with their scale depen- 
dence. The wave function of the 5°-meson, the distribution amplitudes for the light pseudo-scalar 
and vector mesons and the input values of the various mesonic decay constants are also given here. 
Section III contains the calculation of the B® — > PP mesons, making explicit the contributions 
from the so-called emission and annihilation diagrams. Numerical results for the charge-conjugated 
averages of the decay branching ratios, direct CP-asymmetries, the time-dependent CP asymme- 
tries Sf and the observables Hf in the time-dependent decay rates are tabulated in Tables III, 
V and VI, respectively. These tables also contain detailed comparisons of our work with the cor- 
responding numerical results obtained in the QCDF and SCET approaches, as well as with the 
available data. Section IV contains the numerical results for the decays B® — > PV. They are 
presented in Tables VII, VIII and IX for the charge-conjugated averages of the decay branching 
ratios, time-integrated CP-asymmetries, the time-dependent CP asymmetries Sf and the observ- 
ables Hf in the time-dependent decay rates, respectively. We also show the corresponding results 
from the QCDF approach in Tables VII and VIII. Section V is devoted to a study of the decays 
B® — > VV, making explicit the amplitudes for the longitudinal (L), normal (N) and transverse (T) 
polarization components of the vector mesons. Numerical results for the CP-averaged branching 
ratios are presented in Table X, and compared with the corresponding results from the QCDF 



approach, updated recently in Ref. 



181 ]. and available data. Results for the three polarization 



fractions fo, f\\ , f±, the relative strong phases <j)\\ (rad), (f>±(rad) and the CP-asymmetries are 
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displayed in Table XI. Appendix A contains the various functions that enter the factorization for- 
mulae in the pQCD approach. Appendix B gives the analytic formulae for the B® — > PP decays 
used in the numerical calculations, while the details of the formulae for the decays — > PV and 
B® — > VV are relegated to Appendix C and D, respectively. 

II. THE EFFECTIVE HAMILTON! AN AND THE INPUT QUANTITIES 



A. Notations and Conventions 



We specify the weak effective Hamiltonian 



Keff = { V ub V* 
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uq 



C^QUv) + CMQ U M\ - V tb V* ^diriQiifi) 

i=3 

where q = d,s. The functions Qi (i = 1, 10) are the local four-quark operators: 
• current-current (tree) operators 



Ql = (u a 6 / g)y_ A (g/3M Q )y- 



QCD penguin operators 

Q 3 = 



{q a b a ) v . A J2W^ 
g' 

Q 5 = {qab a )v-A^2{q' p q'p 



V-A, 



V+A, 



Ql = (Uab a )v-A(qpUp)v-A, 



Qa = {q/3b a )v-A^{qWp)v-A, 



H.c, 



Qt 



a V-A 



V+A, 



(2) 



(3) 

(4) 
(5) 



electro-weak penguin operators 



Qi = ^{qabajv-A^eq'iq'pq'^v+A, Qs = -(^Mv-A ^ MC^Wa, 



q' 



Q9 = ■z(qab a )v-A^2e q >(q' l3 tfp)v-A, Qw = Aqpb a )v-A ^2 e q'^' a ^)v-A, 



(6) 
(7) 



where a and (3 are the color indices and q' are the active quarks at the scale mj, i.e. q' = 
(u,d, s,c,b). The left handed current is defined as (q^q'^v-A = q'aluiX — lb)q'p an d the right 
^q'^v+A = q' a 7iy(l + ls)q'/3- The combinations a, of Wilson coefficients are 



handed current ' 
defined as usual 



ai = C 2 + d/3, a 3 = C 3 + C4/3, a 5 = C 5 + C 6 /3, a 7 = C 7 + C 8 /3, a 9 = C 9 + do/3, 
a 2 = C 1 + C 2 /3, a 4 = C 4 + C 3 /3, a 6 = C 6 + C 5 /3, a 8 = C 8 + C 7 /3, a 10 = C 10 + C 9 /3. 
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TABLE I: Numerical values of the combinations of Wilson coefficients denned in the text at different 
scales (fi). 



\i (GeV) 


2.5 


2.0 


1.5 


1.0 


a\ 


f.f 


1.1 


1.1 


1.1 


a 2 (xirr 2 ) 


f.f 


-2.8 


-8.7 


-19.4 


a 3 (xl0~ 3 ) 


6.2 


7.5 


9.7 


14.4 


a 4 (xl(T 3 ) 


-32.0 


-35.8 


-41.4 


-51.3 


a 5 (xl0~ 3 ) 


-5.6 


-7.4 


-10.5 


-17.6 


a 6 (xl0~ 3 ) 


-46.8 


-54.9 


-68.2 


-95.6 


a 7 (xl0" 4 ) 


12.6 


12.7 


13.2 


14.0 


a 8 (xlO- 4 ) 


9.6 


10.6 


12.2 


15.7 


a 9 (xlO~ 4 ) 


-84.3 


-85.4 


-87.2 


-91.0 


aio(xfO" 4 ) 


-0.87 


2.3 


7.0 


15.5 



Since we work in the leading order of perturbative QCD (0(a s )), it is consistent to use the leading 
order Wilson coefficients. The scale \i characterizes the typical scale for the hard scattering, and 
for four different values of this scale the results of the Wilson coefficients are listed in Table [H 
It is seen that the coefficient a 2 and the penguin operator Wilson coefficients a 3 , a§ have a large 
dependence on the scale which will give large uncertainties to those channels highly dependent on 
them, such as the color-suppressed QCD penguin dominant processes. This situation is typical of 
leading order estimates; the reduction in the scale-dependence requires the calculation of next-to- 
leading order results which are beyond the theoretical accuracy to which we are working. 

We consider the B° s — > M 2 M 3 decay. In the emission type diagrams with the spectator quark 
s from the initial state B s recombining to form the meson M 3 , we denote the emitted meson as 
M 2 while the recoiling meson is M 3 . The factorization formula is then denoted as Fg s _> M3 . It is 
convenient to define two light-like vectors: n and v. These two vectors satisfy n 2 = v 2 = and 
n-v = 1. The meson M 2 is moving along the direction of n = (1, 0, Ot) and M3 is on v = (0, 1, Ot), 
we use Xi to denote the momentum fraction of the anti-quark in each meson, ki± to denote the 
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FIG. 1: The notation in our calculation 



transverse momentum of the anti-quark: 

Pb = %(M,0 T ), P 2 = ^|(1,0,0 T ), P 3 = ^|(0,1,0 T ), 

M b M b M b 

h = (xi— ^,0,ki_L), k 2 = (x 2 — ^,0, k 2 _L), fc 3 = (0,x 3 — ^,k 3± ), (9) 

which are shown in Figure [U 



B. Wave Functions of the Meson 



In order to calculate the analytic formulas of the decay amplitude, we use the light cone wave 
functions $m,«/3 decomposed in terms of the spin structure. In general, $M,a/3 with Dirac indices 
a, j3 can be decomposed into 16 independent components, l a/ g, 7^, cx^, (7^75) 750/3. If the 
considered meson M is the B s meson^ heavy pseudo-scalar meson, the B s meson light-cone 



matrix element can be decomposed as 



24 



d*ze^ z (0\b a (0)sp(z)\B s (P Bs )) 



(Pb,+M b ,)j 5 



<t>B s {h) - ^ ^ Ba (fci) 



(10) 



Pa 



From the above equation, one can see that there are two Lorentz structures in the B s meson 
distribution amplitudes. They obey the following normalization conditions 



d 4 h 



<pB s (h 



f. 



d 4 h 



<pB s (h 



0. 



11] 



(2vr )4 — V ^ 2%/ g> J ( 27r )4- 

In general, one should consider these two Lorentz structures in the calculations of B s meson 
decays. However, it is found that the contribution of <pB g is numerically small 26j, thus its 
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contribution can be neglected. With this approximation, we only retain the first term in the 
square bracket from the full Lorentz structure in Eq. (fTUl) 



$ Bs = ^=(P Bs + Ms^sXh). (12) 

In the next section, we will see that the hard part is always independent of one of the kf and/or 
ki, if we make the approximations shown in the next section. The B s meson wave function is 
then a function of the variables (or kf) and k± only, 

/dk + 
^ B Xkt,K,ki) . (is) 

Then, the B s meson's wave function in the 6-space can be expressed by 

$b,M) = ^=[p Balb + M Bal ,}<p Bs {x,b), (14) 

where b is the conjugate space coordinate of the transverse momentum k^. 

In this study, we use the model function similar to that of the B meson which is 



4> Bs (x,b) = N Bs x 2 (l -x) 2 exp 



(15) 



with N Bg the normalization factor. In recent years, a lot of studies have been performed for the 
Bj and B^ decays in the pQCD approach Q|. The parameter cjj, = 0.40 GeV has been fixed there 
using the rich experimental data on the B® and B ± mesons. In the SU(3) limit, this parameter 
should be the same in B s decays. Considering a small SU(3) breaking, the s quark momentum 
fraction here should be a little larger than that of the u or d quark in the lighter B mesons, since 
the s quark is heavier than the u or d quark. The shape of the distribution amplitude is shown 
in Fig f5] for uj b = 0.45 GeV, 0.5 GeV, and 0.55 GeV. It is easy to see that the larger oj h gives a 
larger momentum fraction to the s quark. We will use ujf, = 0.50 ± 0.05 GeV in this paper for the 
B s decays. 

C. Distribution Amplitudes of Light Pseudo-scalar Mesons 

The decay constant fp of the pseudo-scalar meson is defined by the matrix element of the axial 
current: 

(0\q llfll5 q 2 \P(P)) = if P P ll . (16) 
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FIG. 2: B s meson distribution amplitudes. The solid-, dashed-, and tiny-dashed- lines correspond to 
u B = 0.45 GeV, 0.5 GeV, and 0.55 GeV. 



TABLE II: Input values of the decay constants of the pseudo-scalar and vector mesons (in MeV) 27j, |28(] 



T 



f. 



f? 
■J u. 



fK* 



Jk* 



f 1 

J A 



fir JK fp Jp Juj J u; J K" J X* 

131 160 209 ± 2 165 ± 9 195 ± 3 145 ± 10 217 ± 5 185 ± 10 231 ± 4 200 ± 10 

The pseudo-scalar decay constants are shown in table [Til taken from the Particle Data 
Group 28|. The vector meson longitudinal decay constants are extracted from the data on 
t~ — > (p~, K*~)v T 28J and the transverse decay constants are taken from QCD sum rules 29i. l30l]. 



The input values given in table [ITJ are very similar to the ones used in 161 ] . 

The light-cone distribution amplitudes are defined by the matrix elements of the non-local light- 
ray operators at the light-like separation with z 2 = 0, and sandwiched between the vacuum and 
the meson state. The two-particle light-cone distribution amplitudes of an outgoing pseudo-scalar 
meson P, up to twist-3 accuracy, are defined by [3jJ: 



(P(P)|g 2 (*)7„76?i(0)|0) 
(P(P)|? 2 (*)75?i(0)|0) 

(P(P)|g 2 (*)<v7 5 gi(0)|0) 



AxP-z . 



-if P P, / dxe tx 
Jo 

—ifprriQ / dxe vxPz (j)^ 
Jo 



2 [X 



-f P mo(PpZ y - P u z ll 
6 jo 



dxe" 



ixP-z la i 



X) 



(17) 
(18) 
(19) 



where we have omitted the Wilson line connecting the two space-time points. <p2{x), <f>£(x) and 
03 (x) have unit normalization. Mp is the mass of the pseudo-scalar meson, m is the chiral 



scale parameter which is defined using the meson mass and the quark masses as mo 



Ml 

m qi +m q2 
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m 



o 



1.4 GeV and =1.9 GeV, x is the momentum fraction associated with the quark q 2 . 



It is easy to observe that the contribution from <p2 {x), independent of the mass, is twist-2. The 
contributions from 4> p {x) and 03 (x), proportional to r = Mq/M^, are twist-3. 
The above definitions can be collected as 



(P(P)\q 2 p(z)q la (0)\0) = ~ [ dxe lxP - z U 5 P^x) +m ol ^ p {x) -m Q a^ lb P,zf 

vo Jo 
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~ J dxe lxP - z [75 F<p A (x) + l 5 m <p p (x) + m o75 (^ ^ - 1)0 T ( 

(20) 



with the redefinitions of the distribution amplitudes: 

<f{x) = <P P (x) = t a (x) = (21) 

and we have performed the integration by parts for the third terms and <p T {x) = ~^</> CT (x). 

In the pQCD approach, the only non-perturbative inputs are the meson decay constants and 
meson light cone distribution amplitudes, and they are both channel independent. The meson 
decay constants are either measured through the leptonic decays of the mesons and the semilep- 
tonic decays of the r-lepton or calculated from the measured ones using broken SU(3) symmetry. 
Therefore there is not much uncertainty in them. The wave functions depend on the factorization 
scale and also the factorization scheme. In principle, they should be determined by experiment. 
Although there is no direct experimental measurement for the moments yet, the non-leptonic B° 
and B ± decays already give much information on them [33I] . Since the pQCD approach gives 
very good results for these decays, especially the direct CP asymmetries in B° — > tx + ti~ and 
B° — > K + ir~ decays 34|, we will use the well constrained light cone distribution amplitudes of 



the mesons in these papers 



311 ] (see 32| for a summary and update of the LCDAs): 



tf{x) = ^x(l -x)[l + 0.44C 2 3/2 (t)], (22) 

<t> P M = ^[l + 0.43C 2 1/2 (t)], (23) 

fi(x) = -^\C\'\t) + 0.55C 3 1/2 (t)], (24) 

<f>£(x) = ¥£x{l -x)[l + 0.17Cf 2 (t) + 0.2C 2 3/2 (t)], (25) 

= ^[l + 0.24C 2 1/2 (t)], (26) 

f K {x) = -i^[cf(f) + 0.35C 3 1/2 (t)], (27) 
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with Gegenbauer polynomials denned as: 

C\ /2 (t) = t, Cf /2 (i) = 3t 

Cl /2 (t) = i(3t 2 - 1), C 2 3/2 (t) = |(5t 2 -l), (28) 
Cl'\t) = ±t(5t 2 - 3) , 

and t — 2x— 1. In the LCDAs for <f>^(x), (p^{x) and 4>^(x), we have dropped the terms proportional 
to C , | //2 ' 3/ ' 2 , and take into account only the first two terms in their expansion, consistently with 
the rest of the LCDAs. These distribution amplitudes are very close to the previous QCD sum 



rule results 31 1. In recent years, there have been continuing updates of the light cone distribution 
amplitudes 35J. However, the changes in the coefficients of the Gegenbauer polynomials do not 
affect our results significantly, as shown in the next section. We also point out that the default 
value of the scale at which the Gegenbauer coefficients are given above is 1 GeV. However, the 
scale of the perturbative calculation in the pQCD approach where these LCDAs enter is typically 
2 GeV. The LCDAs can be scaled up, as the required anomalous dimensions are known. We 
have done this and find that the resulting numerical differences are small. Strictly speaking, these 
differences are part of the NLO corrections. With the current theoretical accuracy, they can be 
absorbed in the uncertainties on the input Gegenbauer coefficients in the numerical calculations. 



As 
Stech 



or the mixing of 77 and r/, we use the quark flavor basis proposed by Feldmann, Kroll and 



361 ] . i.e. these two mesons are made of nn = (uu + dd)/y/2 and ss: 

(29) 





with the matrix, 

(cos 6 — sin < 
I , (30) 
sin cos <; 

where the mixing angle <p = 39.3° ± 1.0°. In principle, this mixing mechanism is equivalent to the 



singlet and octet formalism, as discussed in 36], and the advantage here is that explicitly only 
two decay constants are needed: 

(0\n^n\ Vn (P)) = -|=/„F, 

(0\*fl53\Ti s (P)) = lf s P» • (31) 

We assume that the distribution amplitudes of nn and ss are the same as the distribution 
amplitudes of it, except for the different decay constants and the chiral scale parameters. We 
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use (36] 



f n = (1.07 ± 0.02)/; = 139.1 ± 2.6 MeV, f s = (1.34 ± 0.06)/,, = 174.2 ± 7.8 MeV, (32) 
as the averaged results from the experimental data. The chiral enhancement factors are chosen as 

(33) 



rn 



o — [m^ cos </> + sin — (ny — m^) cos ^ sin i 

^f^n In 



ml 



2m< 



| m 2 , cos 2 + m 2 sin 2 <j) ^— (m 2 , — m 2 ) cos sin <j)\ . (34) 



'/ 



V2f a 



There are gluonic contributions which have been investigated in 371 ]. with the result that these 
parts do not change the numerical results significantly. So, we will not consider this kind of 
contribution in this work. 



D. Distribution Amplitudes of Light Vector Mesons 



We choose the vector meson momentum P with P 2 = M v , which is mainly in the plus direction. 
The polarization vectors e, satisfying P ■ e = 0, include one longitudinal polarization vector 
and two transverse polarization vectors ex- Following a similar procedure as for the pseudo-scalar 
mesons, we can derive the vector meson distribution amplitudes up to twist-3 38 



(V(P,e* L )\q 2f} (z)q la (0)\0) = 1= j\ xe ixPz [M v f L M%)+ A + M v<Pv(x)] aP (35) 



(V(P,e* T )\q 2p (z)q la (0)\0) = ^= J 



(36) 



[M v f T (j) v v (x)+ f T , 
+M v ie PLUpul ^e f T u n p v° '^(x)} a/3 , 

for longitudinal polarization and transverse polarization, respectively. Here x is the momentum 
fraction associated with the q 2 quark. We adopt the convention e 0123 = 1 for the Levi-Civita 
tensor ^ va ^ . 

The twist-2 distribution amplitudes for a longitudinally polarized vector meson can be param- 
eterized as: 



3U x(l-x) \l + alcl /2 (t) 

x(l - x) 1 + a\ K ,Cl /2 {t) + a\ K *Cl /2 {t) 



<Pu(x) 

b K *{x) 



V6' 



y/6 
3f<j> /, 
x/6 



X) 



(37) 
(38) 
(39) 
(40) 
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Here fy is the decay constant of the vector meson with longitudinal polarization, whose values 
are shown in table [III The Gegenbauer moments have been studied extensively in the literatures 



38 



39(, here we adopt the following values from the recent updates 



29 



30 



40J: 



IK* 



0.03 ±0.02, a 



11 -a 11 

2p — u 2w 



0.15 ±0.07, a\ K * = 0.11 ± 0.09, aL = 0.18 ± 0.08, (41) 



and we use the asymptotic form [4l| : 

<f>y(x) = 



3/f 



v t\ 



(j)y(x) 



(42) 



2x/6 ' ' 2v/6 

The twist-2 transversely polarized distribution amplitudes 4 >y h ave a similar form as the Ion 



gitudinally polarized ones in eq. fl37H40|) . with the moments 



2a 



301 



a£ K . = 0.04 ±0.03, at p = a^ 



0.14 ±0.06, a^. =0.10 ±0.08, aL = 0.14 ± 0.07. (43) 



The asymptotic form of the twist-3 distribution amplitudes <p v and (f> v are 

3/y 



<j> v {x) 



ft 2 ), <fy( x ) = 



(44) 



8^ 4^ 
The above choices of vector meson distribution amplitudes can essentially explain the measured 
B — > K*(j), B — > K*p and B — ► pp polarization fractions 
branching ratios. 



41 



42 



43 



44|, together with the right 



E. A Brief Review of the pQCD Approach 

The basic idea of the pQCD approach is that it takes into account the transverse momentum 
of the valence quarks in the hadrons which results in the Sudakov factor in the decay amplitude. 
As an example, taking the first diagram in Fig. [31 the emitted particle M 2 in the decay can be 
factored out (in terms of the appropriate vacuum to M 2 transition matrix element) and the rest 
of the amplitude can be expressed as the convolution of the wave functions 4>b b , 4>m 3 and the hard 
scattering kernel T#, integrated over the longitudinal and the transverse momenta. Thus, 

M oc / dxidx 3 — -^7- -^(j)B(xi,k 1T ,p 1 ,t)T H (x 1 ,x 3 ,ki T , fc 3T ,t)0y(x 3 , k 3T ,p 3 ,t) . (45) 
Jo J (2tt)^ (2tt)^ 

It is convenient to calculate the decay amplitude in coordinate space. Through the Fourier trans- 
formation, the above equation can be expressed as: 

M oc J dx x dx 3 J d%d%(j) B (x 1 ,b 1 ,pi,t)T H (x 1 ,x 3 ,bi,b 3 ,t)(p v (x 3 ,b 3 ,p 3 ,t) . (46) 
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Loop effects can, in principle, be taken into account in the above expression. In general, individual 
higher order diagrams suffer from two types of infrared divergences: soft and collinear. Soft 
divergence arise from the region of a loop momentum where all it's components in the light-cone 
coordinate vanish: 

P = (/+,/-,f r ) = (A,A,A). (47) 

Collinear divergence originates from the gluon momentum region which is parallel to the massless 
quark momentum, 

F = (/+, r , l T ) ~ (m B , A 2 /m B , A). (48) 

In both cases, the loop integration corresponds to J d^l/l 4 ~ log A, so logarithmic divergences are 
generated. It has been shown order by order in perturbation theory that these divergences can be 
separated from the hard perturbative kernel and absorbed into the meson wave functions using 
the eikonal approximation [45] . One also encounters double logarithm divergences when soft and 
collinear momenta overlap. These large double logarithm can be re-summed into the Sudakov 
factor and the explicit form is given in Appendix |A] 

Loop corrections to the weak decay vertex also give rise to double logarithms. For example, the 
first diagram in Fig. [3] gives an amplitude proportional to l/(x\xi). In the threshold region with 
x 3 — > 0, additional collinear divergences are associated with the internal quark. The QCD loop 
corrections to the weak vertex can produce the double logarithm a s In 2 x 3 and the re-summation 
of this type of double logarithms leads to the Sudakov factor St(x 3 ). Similarly, the re-summation 
of a s \n 2 X\ due to loop corrections in the other diagram lead to the Sudakov factor S t {x\). These 
double logarithm can also be factored out from the hard part and grouped into the quark jet 



function 



461 ] . This type of factor decreases faster than any power of x as x — > 0, so it removes 



ras been parameterized in a form which is 



47]. 



the endpoint singularity. For simplicity, this factor 
independent of the decay channels, twist and flavors 

Combining all the elements together, the typical factorization formula in the pQCD approach 
reads as: 

M. oc J dx\dx 3 J d 2 bid 2 b 3 (pB(xi,bi,pi,t) 

xT H (xi,x 3 ,b 1 ,b 3 ,t)4>M 3 (x 3 ,b 3 ,p 3 ,t)St(x 3 )exp[-S B (t) - S 3 (t)] . (49) 

The threshold Sudakov function S t (x 3 ) and the Sudakov exponents Ss(t), S 2 (t) and S 3 (t) are 
given in Appendix A. Again, strictly speaking, the Sudakov improvements result from higher 
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M3 9 M3 E M3 M3 



FIG. 3: The Feynman diagrams for emission contribution, with possible four-quark operator insertions 

order contributions. However, they are included traditionally in the pQCD approach, though in 
most applications, the perturbative function Th(x\, x 3 , bi, 63, t)(j)M 3 {x 3 , b 3 ,p 3 , t) is calculated only 
in the leading order in a s . 

III. CALCULATION OF THE B s -> PP DECAY AMPLITUDES IN THE PQCD AP- 
PROACH 

In the following we will give the general factorization formulae for B s — > PP decays, and the 
pQCD functions can be found in Appendix |A] We will use LL to denote the contribution from 
(V — A)(V — A) operators, LR to denote the contribution from (V — A)(V+A) operators and SP to 
denote the contribution from (S — P)(S + P) operators which result from the Fierz transformation 
of the (V - A)(V + A) operators. 

A. Emission Diagram 

The emission diagrams are depicted in Figure [3j The first two diagrams are called factorizable. 
They will give the B — » M3 decay form factor, if we factor out the corresponding Wilson coefficients 

• (V - A)(V - A) operators: 

fM 2 F^M 3 { a i) = 8irC F M%J M2 J dx x dx 3 J b 1 db 1 b 3 db 3 (j) Bs (x 1 ,b 1 )^ai(t a )E e (t a ) 

'1 + x 3 )<f>i(x 3 ) +r 3 (l - 2x 3 )(0f(x 3 ) +03(x 3 )) h e (xi,x 3 ,bi,b 3 ) 
+2r 3 ^(x 3 )a i (t' a )E e (t' a )h e (x 3 , x h 63, 6 X )}, (50) 

(V - A)(V + A) operators: 

Fb^mM) = -Fb^mM), (51) 
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fM 2 Fj 



(S — P)(S + P) operators: 

= 167rr 2 C F M^jM 2 J dx x dx 3 J b 1 db 1 b 3 db 3 (f) Ba (x 1 ,bi)^a i (t a )E e (t a ) 

x 3 (x 3 ) +r 3 (2 + x 3 )0f(x 3 ) - r 3 x 3 </)3 (x 3 ) /i e (xi, x 3 , 61, 63) 
+2r 3 (f>i{x 3 )a i (t' a )E e {t' a )h e (x 3 ,x 1 ,b 3 ,b 1 )}, (52) 



B^M 3 ( 



with Cf = 4/3 and a,; the corresponding Wilson coefficients for specific channels. In the above 
functions, = m 0i /mB s , where m 0i is the chiral scale parameter. The functions Ei, the definitions 
of the factorization scales tj and the hard functions hi are given in Appendix [A] 

The last two diagrams in Figj3] (c) (d) are the non-factorizable diagrams, whose contributions 

are 

• (V — A)(V — A) operators: 



32ttCfM b /a/6 / dx\dx 2 dx 3 I bidbib 2 db 2 <f)B s {xi,bi)<j) 2 (x 2 ) 



x 



(1 - x 2 )(j> 3 (x 3 ) - r 3 x 3 {(j) 3 {x 3 ) - 03 (x 3 )) ai{t b )E'Jt b ) 



x h n (xi, 1 - x 2 , x 3 , bi, b 2 ) + h n (xi, x 2 , x 3 , b u b 2 ) 



x 



(x 2 + x 3 )0 3 4 (^ 3 ) + r 3 x 3 (0f (x 3 ) + (j> T 3 (x 3 ))] ai(t' b )E' e (t' b )Y (53) 



(V -A)(V + A) operators: 

pi POO 

M B^M 3 ( a i) = 327rC F M| s r 2 /v / 6 / dx x dx 2 dx 3 I bidb 1 b 2 db 2 (j) Bs (x 1 ,b 1 ) 

Jo Jo 



x|/i n (xi, 1 - x 2 ,x 3 ,b 1 ,b 2 ) (1 - x 2 )<f) 3 {x 3 ) (4> 2 {x 2 ) + 2 (x 2 )) 
+^3 (02^(^2) - 02 (^2)) (0f (x 3 ) + 0j (ar 3 )) 
+ (1 - x 2 )r 3 (0f (x 2 ) + 02> 2 )) (03>s) - 0^^)) ] ai(t b )E' e (t b ) 
-h n (x 1 ,x 2 ,x 3 ,bi,b 2 ) x 2 (pf(x 3 )((f) 2 (x 2 ) - 4>l{x 2 )) 
+r 3 x 2 (0f (x 2 ) - 2 P (£ 2 ))(0f (2C3) - 03 (^3)) 

+r 3 x 3 (0f(a; 2 ) + 2 P (x 2 ))(0f (x 3 ) + f 3 {x 3 ))\ ai {t' b )E' e {t' b )^ (54) 
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FIG. 4: The Feynman diagrams for annihilation contribution, with possible four-quark operator inser- 
tions 

• (S — P)(S + P) operators: 

rl poo 

(aj) = 32nC F M^jV6 / dx\dx 2 dx 3 / bidb 1 b 2 db 2 4)B s (xi,bi)4>2(, x 2) 

Jo Jo 

(x 2 -x 3 - 1)03 (x 3 ) + r 3 x 3 (0f (x 3 ) + 03 (x 3 )) 



(55) 



M Bs^M 



x ai (t b )E' e (t b )h n (x u 1 - x 2 , x 3 , b u b 2 ) + a^Q E' e (t' b ) 



x 2 (j)f(x 3 ) + r 3 x 3 (4> 3 (x 3 ) - <f> 3 (x 3 )) h n (x 1 , x 2 , x 3 , h, b 2 



From these formulas we can see that there are cancellations between the two diagrams of Figj3] (c) 
and (d). If the Wilson coefficients are the same, the non-factorizable contributions (proportional 
to the small Xi) are power suppressed compared to the factorizable emission diagram contributions 
in eq. ([501152]). 



B. Annihilation type Diagrams 

We group all the W annihilation and W-exchange, space-like penguin and time-like penguin 
annihilation diagrams together and refer to them as the annihilation type diagrams. In FigjH the 
first two diagrams are the factorizable annihilation diagrams, whose contributions are 
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(V - A)(V - A) operators: 



fB.F t 



LL 



8irC F M^J Ba j dx 2 dx 3 j b 2 db 2 b 3 db 3 {ai(t c )E a (t c ) 

x (x 3 - l)^(x 2 )^(x 3 ) - 4r 2 r 3 <^(x 2 )0f (x 3 ) 
+2r 2 r 3 X30^(x 2 )(^f (x 3 ) - 4%(x 3 )) h a (x 2 , 1 - x 3 , b 2 , b 3 ) 



+ 



x 2 02 (^2)^3 (^3) + 2r 2 r 3 (0f(x 2 ) - ^(x 2 ))^f(x 3 ) 



+2r 2 r 3 x 2 (0f (x 2 ) + <f>l(x 2 ))(j) 3 (x 3 )J a 4 (t' c )E a (t' c )/i a (l - x 3 , x 2 , 63, 6 2 ) } .(56) 

There is a big cancellation between the two factorizable diagrams FigJH(a) and (b), such 
that they are highly power suppressed, which agrees with the long time argument that 
the annihilation contributions are negligible. Especially, if the two final state mesons are 
identical, the formula of eq. (l56l) gives exactly zero. 



(V -A)(V + A) operators: 



F LR (a- 

1 ann\ u i 



(57) 



(S — P)(S + P) operators: 
f B F SP 



I&xCfMbJb. J dx 2 dx 3 J b 2 db 2 b 3 db 3 { 2r 2 0^(x 2 )03(x 3 ) 
+ (1 - x 3 )r 3 0^(x 2 )(0f (x 3 ) + 4> 3 {x 3 )) ai(t c )E a (t c )ha(x 2 , 1 - x 3 , b 2 , b 3 ) 



+ [2r 3 2 4 (x 2 )0f (x 3 ) + r 2 x 2 (<f)%(x 2 ) - ^{x 2 ))(j) 3 {x 3 ) 
xai(t' c )E a (t' c )h a (l - x 3 , x 2 , b 3 , b 2 



(51 



It is interesting to see that the two diagrams FigJH (a) and (b) give constructive contributions 
here. Furthermore, they are not power suppressed as the (V — A)(V — A) operator contribution 
in eq. (1561) (proportional to the small Xi), but proportional to 2r 2 or 2r 3 . This gives the chirally 
enhanced contributions in the annihilation type diagrams. The operator C^-induced space-like 
penguin contributions produce a large strong phase, which is essential to explain the large direct 
CP asymmetry in the B° — > 7r + 7r~ and B° — > K + ir~ decays 34]. 

The last two diagrams in Figj4] are the non-factorizable annihilation diagrams, whose contri- 
butions are 
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(V - A)(V - A) operators: 

rl poo 

M ann(ai) = 32nC F M B jV6 / dx x dx 2 dx 3 \ hdb 2 b 2 db 2 (l) Bs {xxM) 

Jo Jo 



x {h na (xi,x 2 , x 3l 61, 6 2 ) - x 2 (j) 2 (x 2 )(f) 3 (x 3 ) - 4r 2 r 3 0f (a; 2 )0f(a: 3 ) 
+r 2 r 3 (l - x 2 )(^(x 2 ) + 0j(x 2 ))(0f (x 3 ) - $(x 3 )) 
+r 2 r 3 x 3 (0f (x 2 ) - 02 (x 2 ))(0f (x 3 ) + 3 (x 3 )) ai(t d )E' a (t d ) 
+h' na (xi,x 2 , x 3 , 61, 62) (1 - z 3 )0 2 (x 2 )0 3 (ar 3 ) 
+ (1 - x 3 )r 2 r 3 (0f(x 2 ) + <%(x 2 )){<f>£ (x 3 ) - <f%(x 3 )) 
+x 2 r 2 r 3 (^(x 2 ) - 0^(x 2 ))(0f (x 3 ) + ^(xg))!^^)^^) 



(59) 



(V - A)(V + -4) operators: 

dx x dx 2 dx 3 J b 1 db 1 fodb2<f>B.(x 1 ,b 1 ] 



x^h na (x 1 ,x 2 ,x 3 ,b 1 ,b 2 ) r 2( 2 - x 2 )(^f(x 2 ) + (pl(x 2 ))(t) 3 {x 3 ) 

-^(l+Xg)^^)^^)-^^))]^^)^^) 

+/CO1, 2c 2 , a?3, 61, 62) r 2 x 2 (02 (ar 3 ) + 2 (x 2 )) 3 (ar 3 ) 

+r 3 (:r 3 - 1)<#(* 2 ) (0f (* 3 ) - ^(^ 3 ))1 a^K^'A ( 60 ) 



(5 — P)(5 + -P) operators: 

: 32ttCfM b J\/6 / dx\dx 2 dx 3 \ bidbib 2 db 2 (j)B s (xi,bi) 
Jo Jo 



Km 



x |a i (t (i )^(t d )/i na (xi, x 2 , x 3 , 61, 6 2 ) (#3 - ^)^ 2 (x 2 )(j) 3 (x 3 ) 
-4r 2 r 3 0f (x 2 )0f (x 3 ) + r 2 r 3 a? 3 (0f (x 2 ) + (j) 2 (x 2 )) ((j) 3 (x 3 ) - 
+r 2 r 3 (l - x 2 ){^{x 2 ) - <f 2 (x 2 ))(<Pl(x 3 ) + 0^x3)) 
+a i (t' d )E' a (t' d )h' na (x u x 2} x 3 , h, b 2 ) x 2 (f) 2 (x 2 )(j) 3 (x 3 ) 
+x 2 r 2 r 3 (0f (x 2 ) + <j> 2 (x 2 ))((j) 3 (x 3 ) - 4> 3 (x 3 ))) 
+r 2 r 3 (l - x 3 )(0f (x 2 ) - 2 n (x 2 ))(0f (x 3 ) + <%(x 3 )) 



(61) 



are all power suppressed compared to the factorizable emission diagrams. 
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C. Results for B s —>■ PP decays 



First we give the numerical results in the pQCD approach for the form factors at maximal 
recoil. For the form factors, we obtain: 

rpB-t-K _ n 90+O.O5+O.OO T?B^K _ n 90+O.O6+O.OO (nn\ 
r Q — U.ZO_o 04 _o 00, -^o ~~ u - zo -0.05-0.00' Y 3 ^) 

t^Bs^K _ n 9 .+0.05+0.00 

^0 - U -^ 4 -0.04-0.0D ( M ) 

where f B = 0.19 ± 0.02 GeV, u B = 0.40 GeV (for the B ± and B° d mesons) and f Ba = 0.23 ± 0.02 
GeV, ujb s = 0.50 ± 0.05 GeV (for the B° s meson) have been used. They quantify the SU(3)- 
symmetry breaking effects in the form factors in the pQCD approach. The input values for f B 
and f Bs are in agreement with the unquenched lattice results 48j f B = 0.216 ± 0.022 GeV and 
f Ba = 0.259 ± 0.032 GeV, and with the results from the QCD sum rules Q, We also mention 
in passing that a recent calculation of the distribution amplitudes for the pion and kaon with the 



QCD sum rules 5JJ, which includes a new logarithmic divergent term in the twist-3 distribution 



amplitudes, leads to very large SU(3)-breaking effects in the form factors, calculating the form 



factors in the standard pQCD approach 52j. This has also been observed in [27|. Whether this 



feature also emerges in the light-cone sum rule approach is not clear to us. 



53j and drop the (small) 



For the CKM matrix elements, we adopt the updated results from 
errors on V ud , V us , V ts and V tb : 

\V ud \ = 0.974, \V US \= 0.226, \V ub \ = (3.68jft£) x 10~ 3 , 

\Vtd\ = (8.20+°;^) x 10~ 3 , \V U \ = 40.96 x 10~ 3 , \V tb \ = 1.0, (64) 

a = (99±t 4 )°, 7 = (59.0^)°, urg[-V u V£] = 1.0°. 

The CKM factors mostly give an overall factor to the branching ratios. However, the CKM angles 
do give large uncertainties to the branching ratios of some decays and to all the non-zero CP 
asymmetries. We will discuss their effects separately. 

The formulas of the decay amplitudes for the various channels in terms of above topologies 
are given in Appendix [B] The CP- averaged branching ratios of B s — > PP decays are listed in 
Table IIII1 The dominant topologies contributing to these decays are also indicated through the 
symbols T (tree), P (penguin), P EW (electroweak penguins), C (color-suppressed tree), and ann 
(annihilation). The first error in these entries arises from the input hadronic parameters, which is 
dominated by the S s -meson decay constant (taken as f Bs = 0.23 ± 0.02 GeV) and the B s meson 
wave function shape parameter ( taken as u b = 0.50 ± 0.05 GeV). The second error is from the 
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TABLE III: The CP-averaged branching ratios (xlO™ 6 ) of B s — > PP decays obtained in the pQCD 
approach (This work); the errors for these entries correspond to the uncertainties in the input hadronic 
quantities, from the scale-dependence, and the CKM matrix elements, respectively. We have also listed 
the current experimental measurements and upper limits (90% C.L.) wherever available [t|. For com- 
parison, we also cite the theoretical estimates of the branching ratios in the QCD factorization framework 



1J, and in SCET 



201 ]. quoting two estimates in the latter case for some decays. 



Modes 


Class 


QCDF 


SCET 


This work 


Exp. 


B° s - 


>K+TT- 


T 


ln 9+4.5+3.8+0.7+0.8 

iu - z -3.9-3.2-1.2-0.7 


4 Q _|_ 1 0-1-1 5X09 
^±.iy zc i.z zc ±.o zc u.o 


7 R +3.2+0.7+0.5 
' - u -2.3-0.7-0.5 


k n 4- n 7K 4-i 
o.uzcu.iOzci.u 


R° 




C 


n aq+0. 28+0. 22+0. 40+0. 33 
u -* y -0. 24-0. 14-0. 14-0. 17 


0.76 ± 0.26 ± 0.27 zb 0.17 


O 1 r+O. 05+0. 10+0.02 
u - lu -0. 04-0. 05-0. 01 






K+K- 


p 


99 7+3.5+12.7+2.0+24.1 
zz - ' -3.2- 8.4-2.0- 9.1 


_LO . Z ZC U . ( ZC 1 . _L zc u . o 


1 fi +4. 2+7.5+0. 7 
lo ' u -3.2-4.1-0.2 


OA A 4- 1 A 4- A R 

Z^l.4: ZC ±.*± ZC 4:.U 




- l\ S\ 


P 


0A 7 +2.5+13.7+2.6+25.6 
Z*. / _ 2 .4- 9.2-2.9- 9.8 


17.7 zb 6.6 zb 0.5 zb 0.6 


1 K c+5. 0+8. 3+0.0 
id - u -3. 8-4. 7-0.0 






— > it rj 


*EW 


n o7^+ 013 + 030 + 008 + 010 

U.U ( O_ 012-0 .025-0 .010-0 .007 


U.Ul^ ZC U.UU^ ZC U.UUu ZC U.UU^t 

0.016 zb 0.0007 zb 0.005 ± 0.006 


O o^+ - 02 + ' 01 + C1 - 00 
u - uo -0. 02-0. 01-0. 00 


1UUU 


B°- 

13 s 




p nI „ 
iEW 


n -,1+0.02+0.04+0.01+0.01 
u - 11 -0.02-0.04-0.01-0.01 


onfi 4- n nn^ 4- n on9+ 064 

U.UUU ZC U.UUO ZC U.UUZ_o q06 

0.038 zb 0.013 zb 0.016^0 036 


n 1 1 +0.05+0.02+0.00 
u - ii -0. 03-0. 01-0. 00 






-> K°r) 


c 


n oa+0. 19+0. 64+0. 21+0. 16 

— U. ±0 — U.Z f — U.U / — U.Uo 


0.80 ± 0.48 zb 0.29 zb 0.18 
0.59 ± 0.34 zb 0.24 zb 0.15 


n 1 1 +0.05+0.06+0.01 
— u.uo — u.uo — U.Ul 






+ K°r)' 


c 


9 n +0. 3+1. 5+0. 6+1. 5 
z - u -0.3-l.l-0.3-0.6 


4.5 zb 1.5 zb 0.4 ±0.5 
3.9 zb 1.3 zb 0.5 ±0.4 


n 79+0.20+0.28+0.11 
u - ' z -0. 16-0. 17-0. 05 




B° s 


-> 7777 


p 


1 ^ R+l.6+9.9+2.2+13.5 
ld - u -1.5-6.8-2.5- 5.5 


7.1 ±6.4 ±0.2 ±0.8 
6.4 ±6.3 ±0.1 ±0.7 


n +2. 6+4. 7+0.0 
°- u -l. 9-2. 5-0.0 


< 1500 




-»■ 7777' 


p 


n +5. 5+32.4+8. 3+40. 5 
°*- u -5.2-22.4-6.4-16.7 


24.0 ±13.6 ±1.4 ±2.7 
23.8 ±13.2 ±1.6 ±2.9 


91 n +6. 0+10. 0+0.0 
zl - u -4. 6-5. 6-0.0 








p 


a-i 7+4.2+26.3+15.2+36.6 
* ±w -4.0-17.2- 8.5-15.4 


44.3 ±19.7 ±2.3 ±17.1 

49.4 ± 20.6 ± 8.4 ± 16.2 


-ia O+3.2+6.2+0.0 
ll.U_ 2 7 _ 3 9 _o.O 






+ 7T + 7T~ 


ann 


o no/i+ - 003 + 025 + 000 + 163 

U.UZ1_ .003-0.012-0.000-0.021 




n r 7 +0. 16+0. 09+0. 01 
u - ' -0.13-0.10-0.00 


< 1.36 




->7r 7r° 


ann 


n m 9+0.001+0.013+0.000+0.082 
U.UlZ^o. 001-0. 006-0. 000-0. 011 




O 90+0. 08+0. 04+0. 01 
u - zo -0. 07-0. 05-0. 00 


< 210 
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hard scale t, defined in Eqs. (Al) - (A8) in Appendix A, which we vary from 0.75t to 1.25i, 

(5) 

and from A.q CD = 0.25 ± 0.05 GeV. The scale-dependent uncertainty can be reduced only if the 
next-to-leading order contributions in the pQCD approach are known. A part of this perturbative 
improvement coming from the Wilson coefficients in the NLO approximation can be implemented 
already. However, the complete NLO corrections to the hard spectator kernels are still missing. 
The third error is the combined uncertainty in the CKM matrix elements and the angles of the 
unitarity triangle. 

For comparison, we also reproduce verbatim the corresponding numerical results evaluated in 
the framework of the QCD factorization (QCDF) 16| . and the ones obtained using the Soft- 
Collinear- Effective-Theory (SCET) j^. For the decays B® — > PP involving an rj- and/or an rf- 
meson, Ref. [^(J quotes two sets of values, which differ in the input values of the SCET parameters 
specific to the iso-singlet modes, describing gluonic contributions to the B — * rf'> form factors 
and the gluonic parts of the charming penguins. They are not too different from each other 
except for the electroweak-penguin-dominated decay — > ir°r]' '. The errors quoted in the QCDF 
case correspond, respectively, to the assumed variation of the CKM parameters, variation of the 
renormalization scales, quark masses, decay constants, form factors and (whenever applicable) 
the 77-77' mixing angle (collectively called "hadronic 1"), uncertainties in the expansion of the 
light-cone distribution amplitudes (called "hadronic 2"), and estimates of the power corrections. 
The errors shown in the case of the SCET-based results are due to the estimates of the SU(3)- 
breaking, l/m& corrections and the errors on the SCET-parameters. The last column gives the 
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13J, and the upper limits correspond 



current experimental data from the CDF collaboration 
to 90% C.L. 

A number of remarks on the entries in Table III is in order. We note that there is general 
agreement among these methods in the tree and penguin-dominated B® —>■ PP decays, with the 
variations reflecting essentially the differences in the input quantities. This agreement is less 
marked for the color-suppressed and electroweak-penguin dominated decays. For the annihilation 
dominated decays, SCET has no predictions and QCDF has essentially no predictive power as 
indicated by the estimates for the decays B® — > 7r + 7r~ and B® — > 7r°7r°, which vary over more than 
an order of magnitude once the parametric uncertainties are taken into account. First measure- 
ments of the tree-dominated decay B® — > K + tt~ and the penguin-dominated decay B® — > K + K~ 
have been reported by the CDF collaboration, and are in the right ball-park of the predictions 
shown in Table III. A good number of the decays shown in this Table will be measured at the LHC 
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and Super B-factories, which would discriminate among the predictions of the three frameworks. 

The decays B — > tttt, B — > Kir, B s — > Kit and B s — > KK have received a lot of theoretical 
interest, as they can be related by SU(3)-symmetry, a question of considerable interest is the 
amount of SU(3)-breaking in various topologies (diagrams) contributing to these decays. To that 
end, we present in Table IV the magnitude of the decay amplitudes (squared, in units of GeV 2 ) 
involving the distinct topologies: T, V, S, and V^yj for the four decays modes of the B® and 
B® mesons. The two decays in the upper half of this table are related by U-spin symmetry (d — > s) 
(likewise the two decays in the lower half). We note that the assumption of U-spin symmetry 
for the (dominant) tree (T) and penguin (V) amplitudes in the emission diagrams is quite good, 
it is less so in the other topologies, including the contributions from the W-exchange diagrams, 
denoted by 8 for which there are non-zero contributions for the flavor-diagonal states tt + tt~ and 
K + K~ only. The U-spin breaking is large in the electroweak penguin induced amplitudes Vsw-, 
and in the penguin annihilation amplitudes Va relating the decays B d — > K + tt~ and B s — > K + K~ . 
In the SM, however, the amplitudes Vsw are negligibly small. 

TABLE IV: Contributions from the various topologies to the decay amplitudes (squared) for the four 
indicated decays B® — > tt + 7t~; K + 7r~ and B® — > i: + K~; K + K~ . Here, T is the contribution from 
the color favored emission diagrams; V is the penguin contribution from the emission diagrams; £ is 
the contribution from the W-exchange diagrams; Va is the contribution from the penguin annihilation 
amplitudes; and Veyv is the contribution from the electro-weak penguin induced amplitude. See text for 
their definitions. 



mode (GeV 2 ) 


\T\ 2 


\V\ 2 


\£\ 2 


\Va\ 2 




\v £w \ 2 




B d -> 7T+7T- 


0.8 


4.8 x 10" 3 


5.5 x 10" 3 


1.6 x 10' 


-3 


0.6 x 10" 


cs 


B s TT + K~ 


1.0 


5.4 x 10~ 3 





3.3 x 10~ 


-3 


0.8 x 10~ 


-6 


B d -» K + ir- 


1.2 


10.2 x 10~ 3 





2.3 x 10" 


-3 


2.9 x 10- 


-6 


B s K+K- 


1.5 


11.3 x 10~ 3 


3.5 x 10~ 3 


7.3 x 10- 


-3 


0.5 x io- 


■6 



The direct CP asymmetry of B s — > / is defined as 

Adir = BR(B° B -> f) - BR(B» -> f) = \A(B s ^f)\ 2 -\A(B s ^f)\ 2 
CP BR(B° S —> f) + BR(B Q S —> /) \A(B.-+f)\* + \A(B.^f)\* m 1 J 

The numerical results for the direct CP asymmetries in B — » PP decays are shown in Table El 

The first error is from the B s meson wave function parameter uj^ = 0.50^q q5 GeV. As fs s is an 
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overall factor, it drops out in the ratio and hence does not give any uncertainty in the estimates 
of the direct CP asymmetries. The second error is from the hard scale t varying from 0.75£ to 
1.25i (not changing % — 1,2, 3) and Aq^ d = 0.25 ± 0.05 GeV. The third error is again from 
the combined uncertainty in the CKM matrix elements and the angles of the unitarity triangle. 
The observed CP asymmetry Aq P (B® — > K + tt~) by the CDF collaboration is also shown in this 
table, and is found to be in agreement with the pQCD predictions. 

Within (large) theoretical errors, the observed CP asymmetry is also in agreement with the 
SCET estimate but in stark disagreement with the QCDF prediction. This deserves a comment. 
As is well-known, the decay rates and CP-asymmetries in the decays B® — > 7r + 7r - and B® — > K + tt~ 
are related by SU(3) symmetry. The predicted CP asymmetry A^p(B® — »• K + ir~) in the QCDF 
approach shown in Table V is very similar to the SU(3)-related CP asymmetry Aq P (B® — > 7r + 7r~). 
These CP-asymmetries are small in the QCDF approach and both are consistently in disagreement 
with the data, in magnitude and sign. Using SU(3)-symmetry, the CP-asymmetry Aq P (B^ — > 
K + 7i~) in SCET, given in Table V, is identical to the values quoted in Ref. [3] for the CP- 
asymmetry Aq P (B® — > 7r + 7r~). As the dynamical hadronic quantities in the SCET framework 



are fitted using the B-factory data on B — > nix and B — > Kir [20], it is not surprising that this 
phenomenological fit can account for the SU(3)-related observed CP-asymmetry in the B® decays, 
such as Aq P (B® —> K + rc~). The fitted quantities, namely the SCET-specific functions £, (j, differ 
from their corresponding QCDF analogs. But, crucially, also the contributions of the charming 
penguin denoted as |y4 cc | and argA cc in 20], which are power-suppressed in QCDF, are found 



to be large using data and SCET. These differences lead to a different phenomenological profile of 
the B — > hih,2 decays in SCET and QCDF. However, one has to stress that the charming penguin 
contribution, which is not yet shown to factorize in SCET, obviously is not on the same theoretical 
footing as the rest of the decay amplitudes, for which factorization is proven in the heavy quark 



20j essentially reflect the parameterization of the charming 



limit. Thus, the predictions in SCET 
penguins from the available data. 

In general, one has also to take into account the uncertainties caused by the Gegenbauer 
moments. In recent years, the light cone distribution amplitudes have been continually updated 



35j . In order to check the sensitivity to the values on the Gegenbauer moments, we calculate the 



branching ratio for B s — > K + ir with the following values for the twist-2 LCDAs: 

af = 0.06 ±0.03, <% K = 0.25 ±0.15, (66) 
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instead of af = 0.17 and af = 0.2 from eq. (25) and a\ = 0.44 from eq. (22). Using the above 
values and the asymptotic forms for the twist-3 LCDAs, we obtain 

BK(B° S -> K+«-) = (7.2l|°t°;™) x 10" 6 , (67) 

where the errors are to be interpreted as before. The agreement between this result with the 
corresponding number listed in Table III II confirms our expectation that this branching ratio is not 
changed significantly. In our calculation, the parameters in the LCDAs are chosen at fi — 1 GeV 
but as they are scale-dependent quantities, their value at fi = 2 GeV, the typical scale which enters 
the perturbative calculation in pQCD, is required. At \i — 2 GeV, the values for the Gegenbauer 
moments scaled from eq. ( 1661) are: 

af = 0.05 ±0.02, al' K = 0.17 ±0.10. (68) 

Using these values, we obtain: 

BK(B° S -> K+7T-) = (6.9l|°t°£a5) * 10" 6 . (69) 

We remark that the uncertainties caused by the scale dependence of the Gegenbauer moments are 
not large, compared with other uncertainties. 



D. The Observables Sf and Hf in time-dependent decays B®(t) — ► / 



Restricting the final state / to have definite CP-parity, the time- dependent decay width for the 
B s —> f decay is Q : 



T(B»(t) - /) = e- Tt T(B S - /) 



cosh ^— J +H f smh J^ — 



-^cp cos(Amt) - S f sin( Ami) 



(70) 



where Am = m# — > 0, T is the average decay width, and Ar = — Tl is the difference of 
decay widths for the heavier and lighter B® mass eigenstates. The time dependent decay width 
r(Bg(t) — > f) is obtained from the above expression by flipping the signs of the cos(Amt) and 
sin(Amt) terms. In the B s system, we expect a much larger decay width difference (ATJT)b s - This 
is estimated within the standard model to have a value (Ar/r)# s = —0.12 ± 0.05 55J, updated 
recently in 3] to (Ar/r) Bs = -0.147 ± 0.060, while experimentally (AT/T) Bs = -0.33; 



-0.09 
-0.11 



57j, 
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TABLE V: The direct CP asymmetries (in %) in the B s — > PP decays, obtained in the pQCD approach 
(This work); the errors for these entries correspond to the uncertainties in the input hadronic quantities, 
from the scale-dependence, and the CKM matrix elements, respectively. The only measured CP asym- 
metry is also given For comparison, we also cite the theoretical estimates of the CP asymmetries in 



the QCD factorization framework [la ], and in SCET 



Q, 



20] 



Modes 


Class 






T 


Tf,- 




C 


Tf,- 


K + K~ 


P 


Tf,- 


>K°K° 


P 


Tf, 


— > 7r°r/ 


Pew 


Tf, 


-»7rV 


Pew 


«■ 


-► K°7] 


C 



QCDF 



SCET 



This work 



EXP 



-6.7 
41.6 



+2.1+3.1+0.2+15.5 
2.2-2.9-0.4-15.2 

+16.6+14.3+ 7.8+40.9 
12.0-13.3-14.5-51.0 



4.0 



+1.0+2.0+0.5+10.4 
-1.0-2.3-0.5-11.3 



n q+0. 2+0.2+0. 1+0.2 
u - y -0. 2-0.2-0. 1-0.3 



97 q+6. 0+9.6+2. 0+24. 7 
Ll -°-7.1+5. 7-2. 0-27.2 



20 ± 17 ± 19 ±5 

76 ± 26 ± 27 ± 17 
-6±5±6±2 
< 10 



24.1 



+3.9+3.3+2.3 
-3.6-3.0-1.2 



39 ± 15 ± 



cq 4+1.8+ 7.4+2.2 
dy -*-4.0-11.3-3.5 



-23.3 



+0.9+4.9+0.8 



-0.4 



0.2-4.4-1.1 


+0.6+2.2+0.0 
-0.7-2.2-0.0 



90 R+0.0+2. 0+2.8 
zu -°-0.7-2.5-1.2 



46.8 



+18.5+28.6+5.2+34.6 
-13.2-32.2-12.5-45.6 



B s -> K u 77 



P s -> 7777 



5 S -» 7777 



P 



P 



P 



-36.6 



+8.6+6.0+3.8+19.3 
1.2-7.4-2.5-17.3 



1 R +0.5+0.6+0.4+2.2 
" 1 -O-0.4-0.6-0.7-2.2 



n ^+0. 1+0.3+0. 1+0.4 
u -^-0.1-0.3-0. 1-0.3 



2.1 



+0.5+0.4+0.2+1.1 
-0.6-0.4-0.3-1.2 



-56 ± 46 ± 14 ± 6 
61 ±59 ±12 ±8 
-14±7±16±2 
37 ± 8 ± 14 ± 4 
7.9 ± 4.9 ± 2.7 ± 1.5 
-1.1 ±5.0 ±3.9 ± 1.0 
0.04 ±0.14 ±0.39 ±0.43 
2.7 ±0.9 ±0.8 ±7.6 
0.9 ± 0.4 ± 0.6 ± 1.9 
-3.7 ±1.0 ±1.2 ±5.6 



56.4 



+2.9+6.8+3.1 
-3.4-8.0-3.4 



19.9 



+1.6+5.1+1.4 
-1.4-5.0-0.9 



n fi +o. 2+0.6+0.0 

-W.D_g .2-0.5-0.1 



1 q+0. 0+0.1+0.1 
' 1 -°-0.0-0.2-0.1 



1.9 



+0.2+0.3+0.2 
-0.2-0.4-0.1 



"5° + - 



K - 7r°7T 



aim 



aim 



■1.2 



-1.2 



+0.1+1.2+0.1 
-0.4-1.2-0.1 

+0.1+1.2+0.1 
-0.4-1.2-0.1 



so that both Sf and if/, can be extracted from the time dependent decays of B s mesons. The 
definition of the various quantities in the above equation £1X6 cLS follows: 

2/m[A] 2ge[Aj 

5/ -TTprF' / = TTp^' (71) 



27 



with 



A 



2k 



A(B S -> /) 



(72) 



A(B S -> /) ' 

where 77/ is +1(— 1) for a CP-even (CP-odd) final state / and e = axg[-V cb V ts V* s V£\. With the 
convention arg[V c b] = arg[Vy = 0, the parameter can be reduced to e = arg[— Vf S V^]. The results 
of our calculations for the decays A° — > PP are listed in Table [VI] and compared with the ones 



obtained in SCET 



20|. 



E. Specific Tests of the pQCD predictions in P° -> PP and P° -> PP Decays 

In this subsection, we confront the predictions of the pQCD approach to available data in 
the decay modes B® — > PP and B d — > PP, in terms of the ratios of the branching ratios and 
CP-asymmetries. Restricted by the currently available data, we shall confine ourselves to the 
decay modes A° -> K+K~, A° -> A'+vr", A° -> tt+tt- and A° -> A"-tt+. The ratio of the 
branching ratios defined as Ai = BR ( £ ^> K+K ) h as Deen studied at some length in the literature. 
In R\, the numerator is dominated by the penguin amplitude, but the denominator is a mixture 
of tree and penguin amplitudes, and both the numerator and denominator have been measured 
experimentally. It has been argued that the ratio Ri and the two CP asymmetries Scp{B d — ► 
7T+7T-) and C CP (B° d -> tt+tt-) = -A%{B° d -> tt+tt-) of the -> tt+tt^ channel depend, in the 



SU(3) limit, on only two quantities |58j] which can be determined from data and compared with 
the various dynamical models to get a clear picture of two-body non-leptonic decays. However, 
not invoking the SU(3) limit, this system of observables has too many unknowns and a clean test 
of the dynamical models is bogged down in the details of the hadronic input. To enable getting 
cleaner theoretical handles on the underlying dynamics, we calculate the ratio R 2 = ^r^C^+k-) ' 
the ratio A3 and the quantity called A, defined later in this subsection. The last two (i.e. A3 and 



A) have been advocated by Lipkin 60j invoking earlier work by Gronau 59] as precision tests of 
the SM. 

The CDF collaboration has measured the branching ratio of A° — > K + K~ in the form [13I ]: 
U ■ BR(B° -> A+A- 



Using the results 



f d ■ BR(B° d -> Tr+A-) 



0.324 ±0.019 ±0.041. (73) 



f s = (10.4 ± 1.4)%, U = (39.8 ± 1.0)%, (74) 
AA(A° -> tt+tt-) = (5.2 ±0.2) x 10" 6 , (75) 
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TABLE VI: The mixing-induced CP asymmetries St (the first row of each decay channel) and the 
observables Hf (the second row) in the B s — > PP decays calculated in the pQCD approach (This work). 
The errors for these entries correspond to the uncertainties in the input hadronic quantities, from the 
scale-dependence, and the CKM matrix elements, respectively. Estimates of these quantities in SCET 
are also given, quoting two of these for the decays involving an 7/ and/or an r/-meson. 



Modes 


SCET Theory I 


SCET Theory II 


This work 


5° -> K S TT° 


-0.16 ±0.41 ±0.33 ±0.17 
0.80 ±0.27 ±0.25 ±0.11 




n fil +0.08+0.23+0.01 

— o.oo— o.iy— 0.00 

n ro+0. 04+0. 22+0. 03 
u - JZ -0. 06-0. 16-0. 02 


B° -> K~K+ 

S 


0.19 ±0.04 ±0.04 ±0.01 
0.979 ± 0.008 ± 0.007 ± 0.002 




r. 9O+0-03+0.04+0.02 
— 0.03—0.04—0.01 

n qq+0. 01+0. 02+0. 00 

u - ao -0. 01-0. 02-0. 01 


B° s - K°K° 

S 






0.04 
1.00 


B" -> 7T% 
s f 


0.45 ±0.14 ±0.42 ±0.30 
-0.89 ±0.07 ±0.21 ±0.15 


0.38 ±0.20 ±0.42 ±0.37 
-0.92±0.08±0.17±0.15 


n -1 y+0.04+0.10+0.01 
— O.U4— U. 12 — U.U1 

n qq+0. 00+0. 01+0. 00 
u.yy_ 01-0 .02-0 .00 


B° -> ttV 

S 1 


0.45 ±0.14 ±0.42 ±0.30 
-0.89 ±0.07 ±0.21 ±0.15 


0.38 ±0.20 ±0.42 ±0.37 
-0.92±0.08±0.17±0.15 


n 1 7 +0.00+0.07+0.03 

w ■ _L 1 n ni n n n n n 

— U.01 — O.Uo— U.Oo 

n qr+o. 00+0. 01+0. 01 
u - au -o. 00-0. 01-0. 01 


~B° -> K s n 

S O 1 


0.82 ±0.32 ±0.11 ±0.04 
0.07 ±0.56 ±0.17 ±0.05 


0.63 ±0.61 ±0.16 ±0.08 
0.49 ±0.68 ±0.21 ±0.03 


—0 4Q+°- 03 +0- 22 + - 02 

—0.04 — 0.21—0.03 

n 7n +0.04+0.13+0.01 
u - ' u -0. 05-0. 21-0. 01 


B° -> K s n' 

S O 1 


0.38 ±0.08 ±0.10 ±0.04 
-0.92 ±0.04 ±0.04 ±0.02 


0.24 ±0.09 ±0.15 ±0.05 
-0.90 ±0.05 ±0.05 ±0.03 


— KS+ 01 + 06 + 00 
—0.02 — 0.05—0.00 

n 7n+0- 02+0. 06+0. 00 
u - ' u -0. 02-0. 07-0. 00 


B s -> 7777 


-0.026 ± 0.040 ± 0.030 ± 0.014 
0.9965 ± 0.0041 ± 0.0019 ± 0.0015 


-0.077 ± 0.061 ± 0.022 ± 0.026 
0.9970 ± 0.0048 ± 0.0017 ± 0.0021 


n n o+0. 00+0. 01+0. 00 

U-UO_Q QQ_Q .01-0 .00 

1 nn+ - 00 + - 00 + - 00 
i - uu -o. 00-0. 00-0. 00 




0.041 ± 0.004 ± 0.002 ± 0.051 
0.9992 ± 0.0002 ± 0.0001 ± 0.0021 


0.015 ±0.010 ±0.008 ±0.069 
0.9996 ± 0.0003 ± 0.0003 ± 0.0007 


n nA +0.00+0.00+0.00 
u - u ^-o. 00-0. 00-0. 00 

1 nn+ - 00 +°- cl0 + - cl0 
i - uu -o. 00-0. 00-0. 00 




0.049 ± 0.005 ± 0.005 ±0.031 
0.9988 ± 0.0003 ± 0.0002 ± 0.0017 


0.051 ±0.009 ±0.017 ±0.039 
0.9980 ± 0.0007 ± 0.0009 ± 0.0041 


r, n. A +0.00+0.01+0.00 
u - u ^-o. 00-0. 01-0. 00 

1 nn+ - 00 + - 00 + - cl0 
1 - uu -o. 00-0. 00-0. 00 


s° + - 






n -1 a +0.02+0.08+0.09 
u - 1 ^-0. 00-0. 02-0. 05 

n qq+0. 00+0. 00+0. 00 
u - aa -o. 00-0. 01-0. 01 


5° -> 7r°7r° 






n 1/1+0.02+0.08+0.09 
u,1 *-0. 00-0. 02-0. 05 

n qq+0. 00+0. 00+0. 00 
u.yy_ . 00-0. 01-0. 01 



29 



one obtains R\ and R 2 : 



R x = 4.69 ± 0.94, R 2 = 1.24 ± 0.24. (76) 

We have calculated the branching ratios for the related B^ — > 7r + 7r~, tt + K~ decays in the pQCD 
approach, getting: 

BR{W d -> tt+tT) = (5.8^;?+^) * HT 6 , BR(W d - = (ll.e^tl^) x lO^ 6 . (77) 

Taking into account the correlated errors in the numerator and the denominator, we get: 

D _ o qc+0.45+0.99+0.19 d _i i o+0. 11+0. 10+0.01 fjo\ 
I\\ — ^•O0_o.36-0.59-0.15' -"-2 — ± - io -0.11-0.08-0.01- \'°) 

Adding the theoretical errors in quadrature, we get R\ = 2.35lo!7? and R 2 = l.l8~to\l- Hence, in 
the pQCD approach, the ratio R\ is smaller compared to the current data. This originates from 
the fact that BR(B® — > 7r + 7r~) is somewhat larger in the pQCD approach compared to the data, 
and BR(B® — > K + K~) is smaller than the currently measured branching ratio (see Table III). 
Furthermore, the value of R± also depends on the s quark mass through the chiral scale parameter 
tjIq = m 2 K /{m s + m u rf). If we use = (1.9 ± 0.2) GeV, where the errors reflects the uncertainty 
in the quark masses, we find that R\ has an additional uncertainty R\ = 2.35_^49- Adding this 
error with the one quoted above in quadrature yields R\ = 2.35l *g 6 . It would be interesting to 
investigate how the NLO contributions modify the ratio Ri as the tree and penguin amplitudes 



are expected to be renormalized differently including the 0(a%) corrections, as shown in 52| in 



the context of the B d — > PP decays. 

The ratio R 2 , on the other hand, comes out just about right in pQCD, as in this ratio both 
the numerator and denominator are dominated by QCD penguin amplitudes. One also expects 
that this ratio is stable under O(at) corrections, as the denominator in R 2 is itself stable against 



such corrections (see Table III in 52j) and very similar arguments apply to the decay rate for the 
numerator. The stability of the color-allowed QCD penguin amplitudes against one-loop correc- 
tions to the hard spectator scattering is also borne out in the QCD factorization framework 8J]. 
The SU(3)-breaking effects in R 2 are of the same size as the corresponding effects in the form 
factors, typically ± 30%, and this likely is the dominant theoretical uncertainty in this ratio. In 
Fig. [5l we plot the ratio R 2 vs. A^p(K + ir~) for the LO pQCD-based calculations worked out 
by us and compare them with the current data on these observables. The experimental value of 
R 2 has already been given earlier, and we use A^ r P (B d — > K + n~) = —0.093 ± 0.015 Q. We 
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FIG. 5: i?2 vs A^ p (tvK). The region inside the solid (red) box is the pQCD calculation in the LO. The 
experimental results are shown as bands within their ±lcr errors. 

recall that the direct CP asymmetry A^ p (K + n~), calculated by us in the LO pQCD approach, 
a cp(M n + K~) = -OA7t° f 2 -om-om is numerically significantly different than the earlier 
estimates of the same in this approach (for example, the central value of this CP asymmetry in 



LO is quoted as —0.12 in 52J). This mismatch reflects the dependence of A^ p (K + it~) on the 
input quantities, in particular V u b, which have evolved in the meanwhile. Thus, with our input 
values, we find good agreement with data on R2 but not so for Aq P (K + 7t~) , as shown in Fig. [5j 
However, taking into account the O(o^) contributions, this CP-asymmetry is renormalized while 



i?2 remains practically the same. Taking the central values from Table IV of Ref. 52|, one gets a 
K-factor of 0.75 for Aq P (K + tt~) (the central value in NLO is —0.09 with the input values used 
there). Using this K-factor, and our LO calculations, we estimate A^ p (K +/ n-^) = —0.13 ± 0.04 in 
NLO, making it compatible with the current data. This is shown in Fig. [6j It should, however, 



be pointed out for the sake of clarity that the NLO corrections in |52| are not complete, as the 
hard spectator contributions in O(a^) are not all calculated. Hence a residual contribution to 
Acp(K + 7r~) (and other observables) can not be logically excluded. However, our discussion here 
underscores the dominant source of the uncertainty in A < ^ P (K + 7r~), which is of parametric origin. 

In B^ — > K~tt + and —> K + ir~ , the branching ratios are very different from each other due 
to the differing strong and weak phases entering in the tree and penguin amplitudes. However, 



as shown by Gronau 59J, the two relevant products of the CKM matrix elements entering in the 
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FIG. 6: R2 vs Aqp{ttK). The region inside the solid (green) box is obtained by estimating Aqp(ttK) in 
NLO, as discussed in the text. R2 (pQCD) and the experimental results are the same as in the previous 
figure. 



expressions for the direct CP asymmetries in these decays are equal, and, as stressed by Lipkin 60 ] 
subsequently, the final states in these decays are charge conjugates, and the strong interactions 
being charge- conjugation invariant, the direct CP asymmetry in — > K~ir + can be related to 
the well-measured CP asymmetry in the decay B® — > K + ir~ using U-spin symmetry. In this 

3 



symmetry limit, we have 



59 



\A(B S - 7r + /l-)| 2 - \A(B S - 7v-K + )\ 2 = \A(B d - -k + K~)\ 2 - \A(B d - v~ K + )\\ (79) 
A&(B d - n+K-) = -A&(B a - n~K+) ■ BR{Bs - n+K ~] ■ 4^4- (80) 

Following the suggestions in the literature, we can test these equations and search for possible 
new physics effects which would likely violate these relations. To that end, one can define the 
following two parameters (using Eq. (165]) for the definition of CP asymmetry): 



\A(B S -> n+K~)\ 2 - \A(B a -> n~K+ )\ 2 
\A(B d -> 7T-K+)| 2 - \A(B d Ti+K- 
A d j r P (B d ^n+K-) , BR(B a -> n+K') r(B d ) 



(82) 



A% r p (B s ^ir-K+) BR(B d —> tz + K~) t(B s )' 

The standard model predicts R3 — — 1 and A = if we assume ?7-spin symmetry. Since we have 
a detailed dynamical theory to study the SU(3) (and U-spin) symmetry violation, we can check 
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FIG. 7: i?3 vs A: The red (smaller) rectangle is the pQCD estimates worked out in this paper. The 
experimental results with their ±1<t errors are shown as the larger rectangle. 



how good quantitatively this symmetry is in the ratios R3 and A. We find: 

f? — _1 nn+ - 04 +°- 03 + - 09 A — _n nn+ - 03 +°- 03 + - 06 

SI3 — J--UU_ ao4 _ .o4-0.08' ^ ~~ u - uu -0. 03-0. 02-0. 04) \°°) 

The differing values of ujb and ujb s , which enter in the and 5°-meson wave functions, contribute 
dominantly to the first errors, whereas the current uncertainties on the CKM angle 7 and Vu are 
reflected in the third errors given above. The scale-dependent uncertainties leading to the second 
error are relatively small. Adding all the theoretical errors in quadrature, we get R3 = — I.OOIq \q 
and A = -0.00^;^. Thus, we find that these quantities are quite reliably calculable, as anticipated 
on theoretical grounds. On the experimental side, the results for R3 and A are 13]: 



R 3 = -0.84 ±0.42 ±0.15, A = 0.04 ±0.11 ±0.08. (84) 

We conclude that the SM is in good agreement with the data, as can also be seen in Fig. [7] where 
we plot theoretical predictions for R% vs. A and compare them with the current measurements of 
the same. The measurements of these quantities are rather imprecise at present, a situation which 
we hope will greatly improve at the LHC (and Super-B factories). 
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IV. CALCULATION OF B s TO A VECTOR AND A PSEUDO-SCALAR MESON IN 
THE PQCD APPROACH 



A. Decay amplitudes 

In the decays of B s to a vector and a pseudo-scalar meson, only the longitudinal polarization 
of the vector meson can contribute, thus the decay formulas are very similar to those of the 
B s — > PP decays. Since the Lorentz structure of the vector meson wave functions is different from 
the pseudo-scalar case, there are two kinds of emission diagrams in principle. If the emitted meson 
is a vector meson, the decay amplitudes are the same as that of the B s — > PP case, since they 
are both characterized by the B s — > P transition form factors. For the non-factorizable diagrams 
and also those diagrams in which a pseudo-scalar meson is emitted, the distribution amplitudes 
of the pseudo-scalar meson will be replaced by that of the vector meson as follows: 

(f>£(x) -> fa{x), ${x) -> <j>' 2 {x), <pl(x) -> <j>\{x), m 02 -> M 2 , (85) 

where M 2 is the mass of the vector meson. The factorizable emission topology F§^ P does not 
contribute as the vector meson cannot be generated by the scalar or the pseudo-scalar density. 
Furthermore, we must add a minus sign to </>|(x) and 4>\{x) in the annihilation diagram formulas 
F ann and M ann if the vector meson is on the light s quark side. 

B. Numerical results for B s — > PV decays 

The branching ratios of B s — > PV decays are listed in Table IVHl The direct CP asymmetries 
of the B s — > PV decays are given in Table IVIIIl They are compared with the corresponding 
calculations in the QCDF approach [3] , where the sources of the errors in the numerical estimates 
have the same origin as in the discussion of the B — > PP decays. Comparison of the entries in 
Table |VTT] shows that the two approaches give similar results for the tree-dominated decays and the 
QCD penguin- and electroweak-penguin dominated decays (except for those involving an rj or rj' 
in the final states), and they are radically different for the annihilation-dominated decays (the last 
four entries in this table). The branching ratios for B s — > prj and B s — > prj' are of the same order as 
the QCDF results, but there are large differences in B s — *■ u)(4>)r) and B s — > p(<fr)r]'. In QCDF, the 
branching ratios for B s — > <pr], B s — > (prj' and B s — > 07r° (electro-weak penguin dominated process) 
are of the same magnitude, which implies that the color-suppressed QCD-penguins have a reduced 
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value similar to the electroweak penguins. However, in the pQCD approach, the factorization scale 
is low and this dynamically enhances the color-suppressed QCD penguins sizably. 

Predictions for the CP-asymmetries in the two approaches given in Table rvTTTl are, however, all 
quite different. The zeros shown in some cases for the pQCD approach will be lifted on including 
the neglected small subdominant contributions. Also, as opposed to the QCDF approach, we 
are able to calculate the CP asymmetries in the annihilation-dominated topologies. The mixing- 
induced CP asymmetries («S/)b s and the observable {Hf)B 3 , defined in the previous section, are 
shown in Table IIXI Currently, there is no data to confront the estimates given in Tables IVIIj. IVIIII 
and IIXI and again we hope that this will be remedied at the LHC (and Super-B factories). 



V. B s -» VV DECAYS IN PQCD APPROACH 
A. Decay amplitudes 

There are three kinds of polarizations of a vector meson, namely longitudinal (L), normal 
(N) and transverse (T). The amplitudes for a B s meson decay to two vector mesons are also 
characterized by the polarization states of these vector mesons. The amplitudes for the decay 
B s (Pb) Vi(Pii e-L) + VziPzi e 3 At ) can be decomposed as follows: 



a + —^—P^P v n + i—^— ^ vaP P 2a P 30 



M 2 M 3 a a M 2 M 3 
= A L + A N e* 2 (a = T) • e*(a = T) + i^e a ^e* 2a {a)e* w {a)P 2l P 3p , (86) 

where M 2 and M 3 are the masses of the vector mesons V 2 and V 3 , respectively. The definitions of 
the amplitudes A % [i = L, N, T) in terms of the Lorentz-invariant amplitudes a, b and c are 

A L = a e* 2 (L) • e*(L) + j^f4(L) ■ P 3 e*(L) • P 2 , 

A N = a, (87) 

c 



.4 



T — ■ 

r 2 r 3 

The longitudinal polarization amplitudes for the B s — > VV decays can be obtained from those 
in the B s — > PP decays with the following replacement in the distribution amplitudes: 

02(3)0) 02(3) 0), 02(3)0) -> 02(3)0), 02(3)0) ~> 02(3)0), ( 88 ) 

for the emission diagrams, while 

02(3)0) 02(3)0), 0f(3)O) (~)02(3)O), 02(3)0) ^ H02(3)O), (89) 

35 



TABLE VII: The CP-averaged branching ratios (xlO -6 ) of B s — > PV decays obtained in the pQCD 
approach (This work); the errors for these entries correspond to the uncertainties in the input hadronic 
quantities, from the scale-dependence, and the CKM matrix elements, respectively. For comparison, we 
also cite the theoretical estimates of the branching ratios in the QCD factorization framework 



Modes 


Class 


QCDF 


TTiis work 




— 7V*-4- 


T 


o 7+4.6+3.5+0.7+0.8 
°-' -3.7-2.9-1.0-0.7 


7 R +2.9+0.4+0.5 
' - u -2. 2-0.5-0. 3 


B s ~ 


*p K 


T 

A. 


OA c+11. 9+9- 2+1. 8+1.6 
Z4t -°-9.7-7.8-3.0-1.6 


1 7 0+7.7+1.3+1.1 
-0-5.6-1. 6-0. 9 


B s~ 




n 


n 9 r+0. 08+0. 10+0. 32+0. 30 
u - zo -0. 08-0.06-0.14-0. 14 


n n7 +0.02+0.04+0.01 
u - u ' -0.01-0.02-0.01 


B s~ 


7^0 

->p K 


n 


n fi1 +0.33+0.21+1.06+0.56 
u - ol -0. 26-0. 15-0. 38-0.36 


n ns+ - 02 +°- 07 + - 01 

U - UO -0.02-0.03-0.00 


B s ~ 


TV*0 

-> K rj 


c 


n 9r +0. 15+0.49+0. 15+0.57 
u - ZD -0. 13-0.22-0.05-0. 15 


n i 7+0.04+0.10+0.03 
U- 1 '-0.04-0.06-0.01 


B s - 


■» K u r] 


(7 


n 9O+0.04+0.46+0.23+0.29 
U - ZO -0.04-0.24-0.10-0.15 


n nQ+ - 02 +°- 03 + - 01 

u - uy -0.02-0.02-0.01 


B s ~ 




n 


n m +0.20+0.15+0.68+0.40 
u - oi -0. 18-0. 11-0. 23-0. 25 


n i cr+0.05+0.07+0.02 
u ' io -0.04-0.03-0.01 




K K 


p 


4 -1+1.7+1.5+1.0+9.2 
^-l. 5-1.3-0. 9-2.3 


(y n+1.7+1.7+0.7 
°- u -1.5-1.2-0.3 


B s ~> 


7V*-f- TV — 

A A 


p 


c- c+1.3+5.0+0.8+14.2 
d - J -l. 4-2.6-0. 7- 3.6 


A 7+1.1+2.5+0.0 
' -0.8-1.4-0.0 


B s - 


► If U if 


p 


o q+0. 4+1. 5+1. 3+10. 4 
°- y -0. 4-1.4-1. 4- 2.8 


7 o+2. 5+2. 1+0.0 
7-1. 3-0.0 


7?° 

B s - 


7^*0 7>^ 


p 


/I n+0.4+4.6+1.1+13.2 
^• z -0.4-2.2-0.9- 3.2 


A o+0. 7+2.2+0.0 
^•°-0. 7-1.4-0.0 


B s 


— > IT (p 


r EW 


n i 9 +0.03+0.04+0.01+0.02 
u - iz -0. 02-0.04-0. 01-0. 01 


n 1 fi +0.06+0.02+0.00 
u - ±D -0.05-0.02-0.00 


B s 


— ► p°77 


Pew 


n i 7 +0.03+0.07+0.02+0.02 
u - i ' -0.03-0.06-0.02-0.01 


n nR+ - 03 + - 01 + - 00 

u - UD -0.02-0.01-0.00 


7f 
a s 


-s. n<V 


B EW 


n 9 c+0.06+0. 10+0.02+0.02 
u - zo -0. 05-0.08-0.02-0.02 


n i q+0.06+0.02+0.00 
u ' i,:5 -0.04-0.02-0.01 


Tf, 


— > wr/ 


PC 

7 


n nl 9+0.005+0.010+0.028+0.025 
— U.0U4— U.U03 — U.006— O.UOd 


n n/) +0.03+0.05+0.00 
— U.01 — 0.02 — U.0U 


Tf, 


— > w?/ 


P c 


n n9/i+ - 011 + - 028 +°- 077 + - 042 
U.UZ4_ Q 009 _ 006-0 .010-0 .015 


n ^ a +0.18+0.15+0.00 
u - 44l -0.13-0.14-0.01 


Tf, 


-> 077 


p 


n 19+0.02+0.95+0.54+0.32 
u - iz -0.02-0. 14-0. 12-0. 13 


q fi +l. 5+0.8+0.0 
— 1.0 — 0.6 — 0.0 


Tf, 


— > 07/' 


p 


n n^+o.oi+i.10+0.18+0.40 

u - ud -0. 01-0. 17-0.08-0.04 


n -iq+O.06+0. 19+0.00 
u - ly -0.01-0.13-0.00 


Tf, ■ 


-> A"V 


p 


n 97+0.09+0.28+0.09+0.67 
u - z ' -0.08-0.14-0.06-0.18 


n i fi +0.04+0.09+0.02 
u - iD -0.03-0.04-0.01 


rf, 


-+tt u; 


ann 


» 0.0005 


n nn/i+ - 001 + - 000 + 000 
u.uu^_ 00 i_o .001-0 .000 


Tf,- 




ann 


« 0.003 


n 99+0.05+0.04+0.00 
u - zz -0.05-0.06-0.01 


Tf,- 




ann 


« 0.003 


n 9/1 +0.05+0.05+0.00 
u - z *-0.05-0.06-0.01 


Tf,- 


-+7T°p° 


ann 


« 0.003 


n 9 q+0.05+0.05+0.00 
U - ZO -0.05-0.06-0.01 
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TABLE VIII: The direct CP asymmetries (in %) in the B s — > PV decays, obtained in the pQCD approach 
(This work); the errors for these entries correspond to the uncertainties in the input hadronic quantities, 
from the scale-dependence, and the CKM matrix elements, respectively. For comparison, we also cite the 
theoretical estimates of the CP asymmetries in the QCD factorization framework [if]]. 



Modes 


Class 


QCDF 


This work 


B s^> 


7T A 


T 


n R +0.2+1.4+0.1+19.9 
u - u -0. 1-1. 7-0. 1-20.1 


1Q n +2.5+2.7+0.9 
±».U_2.6-3.4-1.4 


B s~ 


>/9 A 


T 


1 r+0. 4+1. 2+0. 2+12.1 
0.4-1.4— 0.3-12.1 


14 9+2.4+2.3+1.2 
J-^.^— 2.2— 1.6— 0.7 


B s~ 


^0 IS*0 
»7T A 


C 


4 c- 7+14.3+13.0+28.4+80.0 
1 -16.0-11.6-28.0-59.7 


47 -1+7.4+35.5+2.9 
• i — 8.7-29.8-7.0 


B s~ 


-*p A 


C 


94 7 +7.1+14.0+22.8+51.3 
' -5.2-12.4-17.7-52.3 


70 4 +6.4+16.2+2.2 
' °- -11.7-47.8-3.9 


B s ~ 


-> A 7/ 


C 


4n n+17.0+24.6+ 7.8+65.9 
<± U-Z - 11.5-30.8- 14.0-96.3 


ki 9+6.2+14.1+2.0 
UJ " -6. 4-12. 4-3. 3 


77° 
B s ~ 


-> A 77 


c 


co fi +16.9+41.4+19.9+44.9 
°"- D -ll. 9-11. 7-13. 9-35. 7 


c-i 1+4.6+15.0+3.2 
,J1 - 1 -6. 6-18. 2-4.1 


77° 
B s ~ 


-> A U) 


c 


ao q+13.6+18.0+30.6+57.7 
^° ,y -13. 4-18. 2-30. 2-49. 3 


t-9 1+3.2+22.7+3.2 
dz,i -0. 0-15. 1-2.0 


77° 
B s 


A A 


p 


9 9+0.6+8.4+5.1+68.6 
z - -0.7-8.0-5.9-71.0 


o fi fi +2.3+2.8+1.3 
ou - u — 2.3— 3.5— 1.2 


B s 


1^*+ is- 
lS. A 


p 


o 1+1.0+3.8+1.6+47.5 
•J- 1 — 1.1-2.6— 1.3-45.0 


rc 0+4.4+8.5 + 5.1 

tJ "- ) - J -4.9— 9.8— 2.5 


B s ~ 


• A A 


p 


i 7+0.4+0.6+0.5+1.4 
x - 1 -0.5-0.5-0.4-0.8 





B s - 


• A A 


p 


n 9+0.0+0.2+0.1+0.2 

u - z -0. 1-0.3-0. 1-0.1 


n 


B s 


— > H (f> 


A h/vV 


97 9+6.1+9.8+2.7+32.0 
z ' ,z — 6.8— 5.6-2.4-37.1 


1 0+0.3+2.1+1.5 

lo -°-0.4— 1.7— 0.7 


77° 
B s 


— >■ p°r] 


± fj w 


97 q+6. 4+9. 1+2. 6+25. 9 
' ' — 6.7-5.7-2.2-28.4 


q 9+I.O+2.8+O.4 
y ' -0.4— 2.7-0.7 


B° 


r '/ 


^ Fj VV 


90 n+6. 1+10. 3+1. 5+24. 8 
°" — 7.5— 6.3—1.8—27.5 


9c 0+1.3+2.8+3.4 
"■—2.0—3.6—1.5 


Tf, 


— ► 


P,C 




Ifi 7+5.8+15.4+0.8 

10./ 09 in 1 -17 

— o.z — iy.i — 1. 1 


Tf, 


— > OJ?/ 


P,C 




7 7+0.4+4.5+9.4 
'•'-0.1-4.2-0.4 


B°, 


-> 017 


p 


,.+2.0+30.1+14.6+36.3 
°-*-2.1-71. 2-44.7-59.7 


1 0+0.0+0.6+0.1 

1,o -0. 1-0.6-0.2 


t, 


-> <f>rf 


p 


R9 9+15.9+132.3+80.8+122.4 
°^ ,z -10.2- 84.2-46.8- 49.9 


7 O+1.5+1.2+0.1 
' -°-0.5-8. 6-0.4 


Tf,- 


-> 


p 


ln n+3.0+4.7+3.7+5.0 
lu - a -2.4-3. 0-4.1-7.5 





Tf 




ann 




n n +0.0+0. 5+0.8 
°- u -5. 2-3. 3-0.4 


Tf- 




ann 




a fi +0.0+2. 9+0.6 
*-°-0.6-3. 5-0.3 


Tf- 




ann 




1 O+0.9+2.8+0.1 
i -°-0.4-3.5-0.2 


Tf- 


^vr p° 


ann 




1 7+0.2+2.8+0.2 
'-0.8-3.6-0.1 
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TABLE IX: The mixing-induced CP asymmetries {Sf)B 3 and (Hf)B s in B s — ► decays obtained in 
the pQCD approach (This work); the errors for these entries correspond to the uncertainties in the input 
hadronic quantities, from the scale-dependence, and the CKM matrix elements, respectively. 



Modes 


Class 


(Sf) B 


V J /Os 




o°Kc 


C 


r, c-y+0.22 + 0.51 + 0.02 

' —0.17—0.39—0.05 


D QR+0- 1 0+0- 4 6+0.04 
U.OD_ 0.13— 0.15— 0.04 


if 


► Kgui 


c 


n fiq+0.09+0.28+0.01 
U.Do n no nil n no 
— u.uy — U. 11 — U.Uz 


n ^7+0.11+0.31+0.02 

U.O i n iq n^fl n no 
— U. lo— U.Oo— u.uz 


if 


71 (p 


A -CjVV 


n n7 +0.01+0.08+0.02 
u - u ' -0.01-0.09-0.03 


n qo+o.oo+o.oi+o.oi 

U.yo_ .00-0.03-0.00 


B°- 

s 


-> (fin 


Pew 


n i c +0.06+0. 14+0.01 
— 0.0b— O.lo— U.01 


n Qo+o.oi+o.oi+o.oo 

— 0.U1 — 0.U3 — 0.U0 


Tf,- 


*fPrf 


Pew 


n i R +o. oo+o. 10+0.04 

O.IU n nri — f) 19 — 05 


n qc+o.oi+o.oi+o.oi 

U.»<J_o 00 — 09 — 09 


Tf,- 


-> a;?? 


P,C 


n n9 +0.01+0.02+0.00 
u - uz -0.03-0.08-0.00 


n qo+o.oi+o.oi+o.oo 

u - yy -0.01-0.06-0.00 


Tf,- 


-> u;?/ 


P,C 


n i 1+0.01+0.04+0.02 
u - 1J --0.00-0.04-0.03 


n QQ+o.oo+o.oo+o.oo 
u - yy -o.oo-o.oo-o.oo 


Tf,- 


-> 07? 


p 


n n^+ - 02 +°- 07 + - 01 

U - UO -0.01-0.20-0.02 


1 nn+ - 00 + - 00 + - 00 
-■■• uu -o.oo-o.oi-o.oo 


Tf,- 


-> 077' 


p 


n nn+ - 00 + - 02 + - 00 

u - uu -0. 00-0. 02-0. 00 


1 nn+ aoo + aoo + - 00 

i - uu -0.00-0.00-0.02 


Tf,- 




p 


-0.72 


-0.69 


Tf,- 


+ 7T° W 


ann 


n Q7+0.00+0.00+0.11 
u - y ' -0.01-0.00-0.02 


n 99+O.O2+O.OO+O.I2 
u ' zz -0.00-0.02-0.29 


Tf,- 


,7rV 


ann 


n i Q+0.00+0.02+0.01 
u - iJ -0. 00-0. 02-0. 02 


n QQ+o.oo+o.oo+o.oo 
u - yy -o.oo-o.oo-o.oo 



for the annihilation diagrams. The factorizable emission topology contribution F B ^ V (i = 
L, N, T) vanish due to the conservation of charge parity. 

The normal and transverse polarization amplitudes for B s — > VV decays are displayed as 
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follows. For the factorizable emission diagrams shown in FigJ3](a) and (b), the formulas are 

fv 2 FBs^v 3 ( a i) = ^C F M%Jv 2 r 2 J dx x dx 3 J b l db l b 3 db 3 (j) Ba (xi,b 1 )[h e {xi,x 3 ,b l M) 



'o Jo 

xE e (t a ) ai {t a M(x 3 ) + 2r 3 <f> v 3 (x 3 ) + r 3 x 3 (<t> v 3 {x 3 ) - <%(x 3 ))] 



+r 3 [<p v 3 {x 3 ) + ^{x 3 )}E e {t' a )a i {t' a )h e {x 3 ,x u b 3 ,b 1 )j, (90) 

fv2 F B^v 3 ( a i) = l ^C F M B J V2 r 2 / dx x dx 3 \ b 1 db 1 b 3 db 3 4>B a {x 1 M)\h e {x 1 ,x 3 ,b 1 M) 

Jo Jo l 



x[<Pi (x 3 ) + 2r 3 ^(x 3 ) - r 3 x 3 (^(x 3 ) - <j>l(x 3 ))]E e (t a )ai(t a ) 
+r 3 [^(x 3 ) + ^(x 3 )]£; e (0«i(0^(x 3 ,x 1 ,6 3 ,& 1 )), (91) 



Fb%M) = Fs^vMi), (92) 
Fb^vM) = °> ( 93 ) 

with i = L,N, T. 

The non-factorizable emission diagrams are shown in FiglH^c) and (d). Their contributions are 
expressed in the formulas given below: 

fl POO 

M Bg l^ V3 (ai) = 327iC F M Ba r 2 /VG / dx\dx 2 dx 3 \ b 1 dbib 2 db 2 (j) B3 (xi-,bi) 

Jo Jo 

x { [x 2 {<j>l{x 2 ) + a 2 (x 2 M(x 3 ) - 2r 3 (x 2 + x 3 ){4>l{x 2 )^ 3 {x 3 ) + (j> a 2 {x 2 )4> a 3 {x 3 ))] 
h n (x u x 2 , x 3 , 6i, b 2 )E' e {t' b )ai(t' b ) 

+ (1 - x 2 )((j) v 2 (x 2 ) + (f) a 2 (x 2 ))(f)J(x 3 )E' e (t b )a i (t b )h n (x u l -x 2 ,x 3 , 61,62)}, (94) 

M^^ V3 (ai) = QAnC F M B r 2 /\fQ J dx\dx 2 dx 3 J bidbib 2 db 2 (fi Bs (xi,bi)^E' e (t b )ai(t b ) 
x [x 2 (<^(x 2 ) + <i) a 2 (x 2 ))4>i"{x 3 ) - 2r 3 {x 2 + x 3 )((f) v 2 (x 2 )^(x 3 ) 
+(p2(x 2 )(p v 3 (x 3 ))] h n {x u x 2 , x 3 , 61, 6 2 ) 

+ (1 - x 2 ) [<t>Z(x 2 ) + ^(a; 2 )]0 3 r (x 3 )^(t fe )a i (t b )/i n (a; 1 , 1 - x 2 , x 3 , 61, 6 2 )}, (95) 

= 647rC F M B jVQ / dx\dx 2 dx 3 \ b\db\b 2 db 2 <\) Bs {x\-,b\) 
Jo Jo 

xr 3 x 3 (t>T{x 2 ){<f) v 3 {x 3 ) - 3 (x 3 )) 
E' e (t' b )a i (t , h )h n {x 1 ,x 2 ,x 3 ,bi,b 2 ) + E' e (t b )ai(t b )h n (xi, 1 - x 2 , x 3 , 61, 6 2 ) L (96) 
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pi /*oo 

327tCfMb s /V6 / dx\dx 2 dx 3 / bidbib 2 db 2 4>B s (xi, bi) r 2 
Jo Jo 

x {^2(02(^2) - (p 2 (x 2 ))(p 3 (x 3 )E' e (t' b )a i (t' b )h n (x 1 , x 2 , x 3 , 61, 6 2 ) 

+ /l n (xi, 1 - X 2 ,X 3 ,b 1 ,b 2 )[(l - X 2 )(4> V 2 (X 2 ) - (j) 2 {x 2 ))4> 3 ( X 3) 

-2r 3 (l - x 2 + rr 3 )(^(x 2 )<^(x 3 ) - fe)^))^^)^)}, (97) 

/•I ^00 

64ttCfM Bs /V6 / dx\dx 2 dx 3 I bidbib 2 db 2 <f) Bs (xi,bi)r 2 
Jo Jo 

x {^2(02(^2) - <p 2 {x 2 ))<pl(x 3 )E' e (t' h )a i (t' b )h n (x 1 , x 2 , x 3 , bi, b 2 ) 

+h n (x 1 , 1 - X 2 ,X 3 ,b 1 ,b 2 )[(l - X 2 )(4> V 2 (X 2 ) - 4> 2 (x 2 ))(pl(x 3 ) 

-2r 3 (l - x 2 + x 3 )(0^(x 2 )0^(x 3 ) - 0^(x 2 )0^(x 3 ))]^(t b )a i (t 6 )}. (98) 

The factorizable annihilation diagrams are shown in FigJH (a) and (b), and the normal polar- 
ization contributions are: 

= -8nC F M B J Bs r 2 r 3 J dx 2 dx 3 J b 2 db 2 b 3 db 3 ^E a (t c )ai(t c )h a (x 2 ,l - x 3 ,b 2 ,b 3 )) 
[(2 - x 3 ) ($(x 2 )$(z 3 ) + $(z 2 )^(z 3 )) + x 3 {<%{x 2 )<f>S(x 3 ) + <t>t{x 2 )</> v 3 (x 3 ))] 

-h a {l - X B , X 2 , 63, h) [(1 + X 2 )(<^(x 2 )<^(x 3 ) + (j>% (^)0 3 (x 3 ) ) 

-(1 - x 2 )(<^(x 2 )<^(x 3 ) + $(z 2 )$(z 3 ))]£a(0*(0}- (99) 

Note, that large cancellations between the two diagrams FigJH (a) and (b) take place, as a result 
of which contributions from these diagrams are suppressed. 
For the transverse polarization, we have 

fB s Fann T ( a i) = ~ Jb s -^ann' T ( a i) 

= -l^C F M B J Bs r 2 r 3 j dx 2 dx 3 j MfcMfcjM^O^ateO + <fe)0 3 Oz 3 )) 
+ (2 - x 3 ){^ 2 {x 2 )(jf 3 {x 3 ) + <p a 2 {x 2 )<Pl{x 3 ))}E a {t c ) ai {t c )h a {x 2 , 1 - x 3 , b 2 , b 3 ) 
+h a (l -x 3 ,x 2 ,b 3 ,b 2 )[(l -x 2 ){<j> v 2 {x 2 )^ 3 {x 3 ) + ct> a 2 {x 2 )ct> a 3 {x 3 )) 
-(1 +x 2 ){<P 2 >{x 2 )<Pl(x 3 ) + r2(x 2 mx 3 ))}E a (t' c ) ai (t' c )y (100) 

We remark that, although the cancellations in this case are not as severe as for the normal 
polarization case, the dominant contributions are still power suppressed by r 2 r 3 , where r 2 ( 3 ) = 
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A/r SP.N , > 



M 



SP,T , 



M 2 (3)/M Bs . For the (S — P)(S + P) operators, we have 



f B .F t 



SP,T, 



^SP,N, 



2fB s Fann" ( a i 

-32-kCfMUb, 



dx 2 dx 3 / b 2 db 2 b 3 db 3 < r 2 ((f) 2 (x 2 ) + (j) 2 (x 2 ))(j) 3 (x 3 ) 
o Jo 1 

xE a (t c )ai(t c )h a (x 2 , 1 - x 3 , h, b 3 ) 

+r 3 ^(x 2 )(^(x 3 ) - </>$(x 3 ))E a (Qai(Qh a (l - x 3 , x 2 , & 3 , & 2 )}. (101) 

Again, for this case, like the 5 S — > PP decays, no cancellations among the contributing diagrams 
or power suppressions are involved. The chiral enhancement here for the transverse and normal 
polarizations are essential for the explanation of the large transverse polarization fraction in the 



penguin dominant B decays, such as B — > K*<p and B — > K*p decays 
factorizable annihilation diagrams shown in FigJD (c) and (d), we have 

M LL,N 
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43j . For the non- 



-MirC F Mg 3 r 2 r 3 /V6 / dx 1 dx 2 dx 3 I b 1 db 2 b 2 db 2 (f)B 3 (x 1 ,b 1 ){(j) 2 , (x 2 )(j) 3 '(x 3 ) 

Jo Jo 

+(j) 2 l (x 2 )(j) 3 l (x 3 )}E ; a (t d )a i (t d )h na (x 1 ,x 2 , x 3 , b u b 2 ), (102) 



I LL ' T (a-) 

l ann l u V 



-M, 



SP,T, 



rl r oo 

-128nC F Mz s r 2 r 3 /\/6 / dx\dx 2 dx 3 / b\db 2 b 2 db 2 ^B s { x U bi)[(i> 2 (x2)(j> 3 (x3) 

Jo Jo 

+<fr 2 (x 2 )(j) v 3 (x 3 )}E' a (t d )a i (t d )h na (x 1 ,x 2 , x 3 , b u b 2 ), (103) 



M LR,T (a- 

IY1 ann \ u i 



= -64:7iC F MgjV6 J dx\dx 2 dx 3 J b 1 db 1 b 2 db 2 4>B s (xi, bi)^h' na (xi, x 2 , x 3 , bi, b 2 ) 

[V2X 2 {4> V 2 {X 2 ) + ^Mix,) - rs(l - X 3 )<^(x 2 )(^(x 3 ) - 03^3))] KVdMQ 

+ [r 2 (2 - x 2 )(ct> v 2 (x 2 ) + r 2 (x2M(x 3 ) - r 3 (l + x 3 )<f 2 (x 2 )(<t>l(x 3 ) - 0f(x 3 ))] 
xE' a (t d )a i (t d )h na (x 1 , x 2 , x 3 , bi, b 2 )\. (104) 

These contributions are all power-suppressed as expected. 



B. Numerical results for B s — > VV decays 



The decay width for B s 



V 2 V 3 is given as 

Pr 



8vrM 



B. 



(i) , 



(105) 



« = 0,||, -1 
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where P c is the momentum of either of the two vector mesons in the final states. The sum is over 
the three transversity amplitudes of the two vector mesons, defined as follows: 

A = -A L , 
A\\ = V2a, 

A ± = r 2 r 3y /2(K? - 1)A T , (106) 

with the ratio k = P2 ■ P3/ '(M2M3). Note that the definitions of Ai(i = 0, ||,-L) are consistent 
with those in 18J, except for an additional minus sign in A Q , so that our definitions of the relative 
strong phases 0j(z =||, _L) (see text below) also differ from the ones in [lij by it. The polarization 
fractions fi(i = 0, ||, _L) are defined as follows: 



* = w + 1^,1, + • (107) 

We first give the numerical results of the form factors at maximal recoil: 



yB- 


>K* n 9^+0.05+0.00 

— U.ZO_Q Q 5 _Q QO, 


A B- 




n on+0.06+0.00 
- U.OU_ 05-O.OI; 


At 


>K* . 


n iq+O.04+0.00 
- U.13_„ 03 _ 00 , 


yB- 


*p n 91 +0.05+0.00 

— U.Z1_Q Q 4 _Q QO, 


A B - 


-*p _ 


n 9 c +0.05+0.00 

u - ZJ -0. 04-0.01; 


A?~ 


>P _ 


n 1 7+0.04+0.00 
u - ± ' -0.03-0.00) 


yB- 


^ n iQ+o.04+0.00 

— u - 1J -0.04-0.00' 


A B- 
^0 




n 90+O.O5+O.OO 
u - zo -0.04-0.01' 


Af~ 


*io _ 


n 1 cr+0.03+0.00 
u ' ld -0.03-0.00- 


V B S 


->K* n 91 +0.04+0.00 

— U.Z1_Q Q3_Q Q1 , 


A B S 




n 9 r +0.05+0.00 
— U.ZO_ 05 _ Q1 , 


Af- 




n R +0.03+0.00 
— U.1U_ 03 _o 01, 


V B S 


^<t> n 9^+0.05+0.00 

— U.ZO_Q 04 _Q Q1 , 


aB 3 


^ = 


n on+O OS+O 00 
- U.OU_ 05 _o 01, 




^ = 


n 1Q +0.03+0.00 

- u.±y_ 03 _o oi- 



(108) 



where the first error in the above entries is due to the input hadronic parameters fB(B s ), &b an d 
the second one is from the hard scale and Aqcd- The entries in the first three lines involve the 
decays of the and mesons to two vector mesons, and as such are not required for the 
B° s — > VV decays being worked out in this paper. We list them here to see the SU(3)-breaking 
effects in the form factors. Within errors, these form factors are in agreement with the ones in 



the light cone QCD sum rules 6JJ and with the slightly different estimates of the same in 



The branching ratios of B s — > VV decays are listed in Table Kl They are compared with the 



corresponding results in the QCDF approach 181 ] . A comparison shows that the tree-, penguin- 



and electroweak dominated decays are comparable in the two approaches, with the numerical 
differences reflecting the input parameters. The color- suppressed decays, but more markedly the 
annihilation-dominated decays, differ in these approaches, a feature which is well appreciated in 
the literature. The experimental upper bounds on some of the decays are also listed. Except for 
the decay B s — > 00, where an experimental measurement may be just around the corner, all other 
decay modes remain essentially unexplored. 
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TABLE X: The CP-averaged branching ratios in B s — > VV decays (xlO 6 ) obtained in the pQCD 
approach (This work); the errors for these entries correspond to the uncertainties in the input hadronic 



quantities, from the scale-dependence, and the CKM matrix elements, respective 
also cite the updated theoretical estimates in the QCD factorization framework 
upper limits (at 90% C.L.) are from the Particle Data Group [28 1. 



y. For comparison, we 
181 ] . The experimental 



Channel 


Class 


QCDF [181 


This work 


Exp [281 


B s - 


■+ p u K* [) 


H 


i c+l.0+3.1 
i -°-0.5-1.5 


n Qo+0.09+0.14+0.00 
U - OJ -0.07-0.09-0.01 


<r 767 


B s - 


-> iL>K* l) 


C 


i 9 +0.7+2.3 
1,z -0.3-l.l 


n o-i +0.10+0.12+0.07 
u - oi -0.07-0.06-0.02 




B s — 


> p K ~ 


T 

J. 


9 c 9 +1.5+4.7 
3 -1.7-3.1 


9n q+8. 2+1.4+1. 2 
zu - y -6. 2-1.4-1.1 




B s - 


K*~K*+ 


P 


q 1+2.5+10.2 
" — 2.2— 5.9 


c 7+1.5+3.4+0.5 
1 —1.2 — 1.4—0.2 




B s ^ 


K*°K*° 


P 


q 1+0.5+11.3 
J - i -0.4-6.8 


7 q+1. 9+3. 8+0.0 
' -°-l. 5-2. 2-0.0 


< 1681 


B s - 


(pK*° 


P 


n .+0.1+0.5 
u -*-0.1-0.3 


n rk+0.16+0.27+0.10 
u - u,j -0.13-0.18-0.04 


< 1013 


B a 




P 


91 g+1.1+30.4 
zl -°-l.l-17.0 


or o+8. 3+16. 7+0.0 
J0 -°-6. 9-10. 2-0.0 


14 ±8 


B s - 


p + p~ 


ann 


n oa +0.03+0.60 
u - M -0.03-0.38 


i n +0. 2+0.3+0.0 
1 - u -0.2-0.2-0.0 




B s 


-^p°p° 


ann 


n i 7+0.01+0.30 
u - 1 ' -0.01-0.19 


n ci +0.12+0.17+0.01 
U - Ji -0.11-0.10-0.01 


< 320 


B s 


^ p°u, 


ann 


< 0.01 


n nn7+ - 002 + - 001 + - 000 
U.UU 1 _o.001-0.001-0.000 




B s 


— > bjijj 


ann 


n i i +0.01+0.20 
u - ii -0.01-0.12 


n oQ+0.09+0.13+0.01 
u - oy -0.08-0.07-0.00 




B s 




Pew 


n ^n+0.12+0.25 
u -* u -0. 10-0.04 


n 90+0.09+0.03+0.00 

u ' zo -0.07-0.01-0.01 


< 617 


B s 


— > 4>uj 


P 


n 1 n +o. 05+0.48 

u - iu -0.03-0.12 


n ic+o.09+0.10+0.01 

u - lo -0.05-0.04-0.00 





The results for the longitudinal polarization fraction fo, parallel polarization fraction f\\ and per- 
pendicular polarization fraction f±, their relative phases <p\\ = Arg(A\\/A ) and = Arg(y4^/y4 ) 
(both in radians), and the direct CP-asymmetry in the B® — > VV decays are displayed in Ta- 
ble IXII In calculating the CP-asymmetries, we have used the definition given earlier in Eq. 
Again, there are no data available to confront the entries in Table IXII 

Two remarks on the B s — > <p<f) decays presented above are in order. 



First, B s — > <fxp is a b — > s penguin dominated process. The CDF collaboration [62| has 
reported the decay branching fraction of this channel as (14 ± 8) x 10 -6 , but a thorough angular 
analysis is still lacking. This channel is very similar to the decay B — > (j)K*, which is well 
measured in the experiment. In the B — > <pK* decays, the data show that the fraction of the 
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left-handed polarization reaches about 50%. This result is quite different from the expectation 
in the factorization assumption that the longitudinal polarization should dominate due to the 
quark helicity analysis. There exist lots of theoretical attempts to solve this contradiction. To 
distinguish, which one is the most appropriate scheme, we must investigate the decay B s — > <fxp 
and other similar channels. 

In the pQCD approach, the weak annihilation diagram induced by the operator 0$ can enhance 
the transverse polarization sizably. In Ref. 63j, the longitudinal fraction of B — > (pK* is about 75 
%, which is much smaller than that obtained based on the factorization assumption, but still larger 



than the data. Moreover, the branching fraction is overestimated. Li [64j has suggested a strategy 
to solve these two problems together, by invoking a smaller value for the form factor A . A smaller 
Aq is also consistent with the predictions in other approaches, such as the covariant light front 
quark model 65J. In the pQCD approach, a smaller A Q requires smaller Gegenbauer moments of 
the longitudinal polarized distribution amplitudes of the recoiling vector meson. Following this 



lead, and performing a calculation using the asymptotic distribution amp 



itude for K*, one indeed 



finds that the longitudinal polarization fraction can be reduced; Ref. 64| finds that this fraction 
can be reduced to 59%. 

Based on the above discussion, we also analyzed the B s — > VV decays adopting the asymptotic 
forms of the twist-3 distribution amplitudes in the pQCD approach, while keeping the leading- 
twist distribution amplitudes up to the second Gegenbauer moments, as done for the pseudo-scalar 
meson case discussed earlier. In addition, we also test the sensitivity of the branching ratio of 
the B s — ► <fi(f) decay on the Gengenbauer moments of the twist-2 distribution amplitudes. The 
result of the branching ratio for this channel with asymptotic twist-2 distribution amplitudes is 
22.2 x 10~ 6 , which is much smaller than that for the case with higher Gegenbauer moments and 
is closer to the experiment. However, we find that the polarization fractions do not change too 
much using the asymptotic forms of distribution amplitudes. The reason is that the contributions 
from the annihilation topology, which enhance the transverse polarization fraction, also decrease 
with a decreased value of the form factor A . This is different from the B — > <pK* decay. 

Second, it is to be noted that the longitudinal polarization of B — > <pK* and B s — > (fxfi can 
be related to each other in the SU(3) limit. The discrepancy between them represents SU(3) 
symmetry breaking effects, which can be reflected in the following aspects in the pQCD approach 
specifically: In the first place, the shape parameter of B and B s wavefunction in the initial 
state can give rise to differences in these two U-spin related process. Then, the longitudinal and 
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TABLE XI: The CP-averaged polarization fractions, relative phases and direct CP asymmetries in the 
B s — ► VV decays obtained in the pQCD approach; the errors for these entries correspond to the un- 
certainties in the input hadronic quantities, from the scale-dependence, and the CKM matrix elements, 
respectively. 



Channel 


/o(%) 




/||(%) 




/x(%) 






</>H (rad) 




<fi±(rad) 


CP(%) 






<r r +0.4+6.9+0.6 

<±J. J_Q.3_ 4,3_ 0.9 


27.6 


fO. 1+2.1- 
-0.2-3.4 


4-0.4 
-0.3 


0R q+0.2+2.3+0.4 
zu - J -0. 3-3. 5-0. 3 


2 


7 


f 0.2+0.2+0.3 

-0.3-0.3-0.1 


2 


0+0.3+0.3+0.4 

°— 0.2-0.2— 0.1 


R1 0+3.2+17.1+4.4 

«J-o_ 4.7-22.8-2.3 




-> loK*° 


ro 9+0.3+3.5+2.3 
JO - z -0.2-2.9-1.3 


23.6 


fO. 2+1.5- 
-0.1-1.7 


t-0.4 
-1.0 


90 -1 +0.1 + 1.4+1.0 
ZO - L —0. 2-1.7-0. 6 


1 


4 


4-0. 1+0.2+0.0 
-0.1-0.2-0.1 


1 


a +0.1+0. 2+0.0 
*-0. 1-0.2-0.0 


R0 1 +4.8+19.7+5.5 
uz - 1 -3.9-12.6-1.9 


B s - 


> p-K*+ 


Q o 7+O.I+O.2+O.O 
y0 - ' -0.2-0.3-0.2 


3.4j 


-0.1+0.24 
0.0-0.1- 


0.1 
0.0 


9 n+0. 1+0. 1+0.1 
z - J -0. 1-0. 1-0.0 


3 





4-0. 1+0. 1+0.1 
-0.0-0.0-0.0 


3 


1 +0.0+0.0+0.0 
1 -0. 1-0. 1-0.1 


9 + I.O+I.2+O.4 
°- z -l. 2-1. 7-1.1 


B s ^ 


K*-K*+ 


ao o+5. 1+2. 1+3.7 
<±C.O_ 4 .0-2.3-1 .5 


30.1 


f 2. 1+0.9+0.8 
-2.7-1.0-1.9 


0(\ 1 +1-8+1.4+0.7 
zu - i -2.4-1.0-1.8 


1 


7 


fO. 2+0. 1+0.1 
-0.2-0.1-0.0 


1 


7+0.2+0.1+0.1 
' -0.2-0.1-0.0 


q q+0.4+3.3+0.3 
J -°-0.7-3.6-0.2 


B s -» 




/IQ 7+5.7+0.6+0.0 
' -4.8-3.8-0.0 


26.8 


4-2.6+2.1- 
-3.0-0.3 


t-o.o 

-0.0 


90 r+2. 2+1. 7+0.0 
zo - d -2. 7-0. 3-0.0 


1 


4 


f 0.2+0.0+0.0 
-0.1-0.0-0.0 


1 


a +0.1+0.0+0.0 
^-0.1-0.0-0.0 





B s - 


-4 (f>K*° 


71 9+3.2+2.7+0.0 
' J -" £ '-3. 0-3. 7-0.0 


15.5 


f 1.6+2. 1- 
-1.7-1.5 


fO.O 
-0.0 


1 q+1. 4+1. 7+0.0 
xo -°-l. 5-1.3-0.0 


1 


4 


+0. 1+0.0+0.0 
-0.1-0.1-0.0 


1 


A +0.1+0.0+0.0 
*-0. 1-0. 1-0.0 





B a 




R1 q+3. 6+2. 5+0.0 
ui - J -3. 2-3. 3-0.0 


20.7 


f 1.7+1.8- 
-2.0-1.4 


fO.O 
-0.0 


-, 7 4 +1.5+1.5+0.0 
1 ' -^-1.6-1.1-0.0 


1 


3 


fO. 2+0. 1+0.0 
-0.1-0.0-0.0 


1 


o+0. 2+0. 0+0.0 
J -o. 1-0. 0-0.0 





B s - 




- 100 




~ 




~ 


4 


3 


4-0. 1+0. 1+0.0 
-0.0-0.1-0.0 


4 


7+0.0+0.3+0.1 
' -0.0-0.5-0.0 


9 1 +0.2+1.7+0.1 
z - i -0. 1-1. 3-0.1 


B s 


-^p a p° 


- 100 




~ 




~ 


4.3 


4-0. 1+0. 1+0.0 
-0.0-0.1-0.0 


4 


7+0.0+0.3+0.1 
' -0.0-0.5-0.0 


9 -1 +0.2+1.7+0.1 
z - i -0. 1-1. 3-0.1 


B s 




- 100 




~ 




~ 


4.5 


f 0.0+0.0+0.0 
-0.0-0.0-0.0 


3 


9+0.0+0.1+0.0 
z -0.1-0.0-0.2 


a n+0. 7+2. 7+1.0 
°- u -0.5-3.9-0.4 


B s 


— > COLO 


- 100 




~ 




~ 


4.3 


fO. 0+0. 1+0.0 
-0.0-0.1-0.0 


4 


7+0.0+0.2+0.0 
'-0.0-0.5-0.0 


9 n +o. 1+1. 7+0.1 
z - u -o. 1-1. 3-0.1 


B s 


^<f>p° 


07 n+0. 2+0.9+0.9 
°' '"-0.2-0.3-0.4 


c o+0. 1+0.2+0.3 
°-°-0. 1-0.5-0.4 


R 9+O.I+O.I+O.I 
u - z -0.1-0.5-0.4 


3.5 


f 0.0+0.0+0.0 
-0.0-0.1-0.1 


3 


r+0. + 0. 1+0.1 

d -0. 0-0. 0-0.0 


-in 1 +0.9+1.6+1.3 

iu - 1 -0.9-1.8-0.5 


B s 


— > <j)UJ 


AA n+0.0+5.4+0.9 
44 - J -7.5-6.1-0.4 


28.5 


4-3.8+3.1- 
-0.0-2.8 


4-0.1 
-0.5 


97 9+3.7+3.0+0.2 
z 1 - z -0.0-2.6-0.4 


3.0 


+0. 1+0.2+0.0 
-0.1-0.2-0.0 


3 


n +0. 1+0.2+0.1 
u -0. 0-0. 2-0.0 


q ^+0.6+2.4+0.6 
°-°-0. 6-2.4-0. 2 



transverse decay constants as well as the Gegenbauer moments in the two vector mesons in the 
final state can also contribute to the SU(3) symmetry breaking effects significantly. In addition, 
these two decay processes also differ due to the absence of the time-like penguin annihilation in 
the B — * <pK* decay. 

VI. SUMMARY 

With the LHC era almost upon us, where apart from the decisive searches for the Higgs boson(s) 
and supersymmetry (or, alternatives thereof), also dedicated studies of the .B^-meson physics (as 
well as that of the heavier 6-hadrons) will be carried out. First results on the decay characteristics 
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of the £?i?-mesons are already available from the Tevatron, in particular the B° s - .B°-mixing induced 
mass difference AM S , the branching ratios for the decays B® — > K + ir~, B® — > K + K~ , B® — »• 00, 
and the first direct CP asymmetry Acp(i?° — > JiT + 7r~). These measurements owe themselves to 
the two experiments CDF and DO; they are impressive and prove that cutting-edge B physics can 
also be carried out by general purpose hadron collider experiments. At the LHC, we will have a 
dedicated 6-physics experiment, LHCb, but also the two main general purpose detectors ATLAS 
and CMS will be able to contribute handsomely to the ongoing research in 6-physics. Conceivably, 
also at a Super-B factory with dedicated running at the T(55')-resonance, the current experimental 
databank on the i? s -meson will be greatly enlarged. 

Anticipating these developments, we have presented in this work the results of a comprehensive 
study of the decays B® — »• /11/12, where hi and hi are light (i.e., charmless) pseudo-scalar and 
vector mesons. This study has been carried out in the context of the pQCD approach, taking 
into account the most recent information on the CKM matrix elements and weak phases and 
updating the input hadronic parameters. The decay amplitudes contain the relevant emission and 
annihilation diagrams of the tree and penguin nature. Explicit formulae for all these amplitudes, 
including the various pQCD-related functions, are provided in the Appendices. Numerical results 
for the charge-conjugation averaged branching ratios, the direct CP asymmetries and the mixing- 
induced CP-asymmetries Sf for the CP-eigenstates /, and the time-dependent observable Hf are 
presented in the form of tables. In addition, for the B® — > VV decays, we also calculate the charge- 
conjugation- averaged transversity amplitudes, their relative strong phases and magnitudes. The 
results for the altogether 49 5° — > h±h2 decays include also theoretical errors coming from the 
input hadronic parameters, variation of the hard scattering scale together with the uncertainty on 
Aqcd; and the combined uncertainty in the CKM matrix elements and the angles of the unitarity 
triangle. For the last mentioned error, we use the recently updated results from the CKMfitter 531 ] . 

Our results are compared with the available data on the decays B® — > h\h 2 from the Tevatron, 
and some selected B® — > h\h% decays from the B-factory experiments, updating the theoretical 
calculations in the pQCD approach with our input parameters. In particular, we revisited the 
well-measured direct CP asymmetry A d ^[B^ d — > K + ir~), working out the parametric sensitivity of 
this observable, and argued that data can be accommodated in the pQCD approach by invoking 
0(a 2 s ) contributions, which we estimated from the existing literature. The successful predictions 
of this direct CP-asymmetry in the pQCD approach is set forth with the first observed CP- 
asymmetry in the decay B® — > K + ir~ . This CP-asymmetry is experimentally large and is in good 



46 



agreement with our numerical results, within the stated errors. We have also analyzed a number 
of ratios of branching ratio involving the B° s — > PP and — > PP decays. They include the ratios 
Ri = BR( - B ^~' K K ) R 2 = BR ( B s-^ K K ) an( j |- wo more called R-x and A, defined in Eq. (72) 
and (73), respectively, which involve the decays B® — > K~ir + and B® — > fT + 7r~ and their charge 
conjugates. Except for the ratio R\, for which the pQCD calculations presented here are on the 
lower side of the data (within large errors), the others are well accounted for. We also compare 
our results with the corresponding ones in the QCDF and SCET approaches. What concerns the 
branching ratios, we find reasonable agreement for some topological amplitudes, but also major 
disagreement, in particular for those decays which are dominated by the color-suppressed tree and 
annihilation amplitudes. There are striking differences in the CP asymmetries, in particular with 
the QCDF-based estimates, which will be precisely tested in the future. 

Experimental precision on B s decays will improve enormously in the coming years, thanks to 
dedicated experiments at the Tevatron and LHC Many of the decay rates and CP asymmetries 
worked out here will be put to experimental scrutiny. They can be combined with the B-factory 
data on the corresponding B — > PP, PV, VV decays, eliminating some of the large parametric 
uncertainties in theoretically well-motivated ratios to test the SM precisely in exclusive hadronic 
decays. In principle, theoretical predictions presented here can be systematically improved by 
including higher order perturbative corrections (in a s ) and sub-leading power corrections in 1/rrib- 
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APPENDIX A: PQCD FUNCTIONS 



In this section, we group the functions which appear in the factorization formulae. 
The hard scales are chosen as 

t a = max{ v / iiM Bs ,l/& 1 ,l/6 3 }, (Al) 
t' a = max{ v ^rM Bs ,l/6 1) l/6 3 }, (A2) 



t b = max{,/x^M Bs , V |1 - x x - x 2 \x 3 M Ba , l/b u l/b 2 }, (A3) 



t' b = m&x{ y /x^M Bs ,y / \x 1 - x 2 \x 3 M Bs ,l/bi,l/b 2 }, (A4) 

t c = max{ VI -x 3 M Ba , 1/hA/h}, (A5) 

t' c = max{V^M Bs ,l/&2,l/& 3 }, (A6) 

t d = max{v/x 2 (l - x 3 )M Ba , ^/l - (l - Xl - x 2 )x 3 M Bs , l/b x , l/6 2 }, (A7) 

^ = max{v/x 2 (l - x 3 )M Bs , vl^i - x 2 |(l - x 3 )M Ba , l/b u l/6 2 }. (A8) 

The functions h in the decay amplitudes consist of two parts: one is the jet function St(xi) 



derived by the threshold re-summation 4a], the other is the propagator of virtual quark and gluon. 
They are defined by 

h e {x u x^b u h) = - h)Io{Vx 3 M B .b 3 )K ( > /£ 3 MB.b 1 ) (A9) 
+9(b 3 - b 1 )I (^ 3 M B b 1 )K (Vx3 M B 3 b 3 )] K (y/xMM B ,bi)St(x 3 ), 
h n (x 1 ,x 2 , x 3 , h, b 2 ) = [6(b 2 - bi)K (^/x^M Ba b 2 )I (^/x^x 3 'M Bs b 1 ) 

+0(&i - b 2 )K (^/xTx^M Bs b 1 )I (^/x^M Bs b 2 )} 
f H^ 1 \^(x 2 -x 1 )x 3 M Bs b 2 ), x x - x 2 < 



K (y/(x 1 - x 2 )x 3 M Ba b 2 ), xi - x 2 > 



(A10) 
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h a (x 2 ,x 3 ,b 2 ,b 3 



h na (x 1 ,X2,x 3 ,b 1 ,b 2 ) = — 



y) 2 ^(x 3 ) [9(b 2 - h)H^\^M B b 2 )J (^/x- 3 M Bs b 3 ) 
+6{b 3 - b 2 )H^\^ 3 M B ^)U^MBM)}H^\V^M B3 b 2 l^l) 
e(h - b 2 )H^\^x 2 {l - x 3 )M Ba b 1 )J (y/x 2 {l - x 3 )M B b 2 ) 



+8(b 2 - b^H^^x^l - x 3 )M B b 2 )JoUx 2 {l - x 3 )M Bs h] 



Ka( X U X 2,X 3 ,bx,b2) 



y 



x K ( a/1 - (1 - xi - x 2 )x 3 M Bs 61 ) , 

6(h - b 2 )H^(y/ X2 (l - X Z )M Ba b 1 )J Q {y/x 2 {\ ~ X 3 )M B b 2 

+6(b 2 - 6 1 )^ 1) ( v /^(l-x 3 )M Bs 6 2 



(A12) 



J 2) J Q 



x 2 {\ - x*)MbM 



X 



(A13) 



JH£ } {^(x 2 - xi)(l - x 3 )M Bs b 1 ), Xl -x 2 <0 
Ko{\/ (xi - x 2 )(l - x 3 )M Bs 6i), xi - x 2 > 

where h£\z) = J (z) +iY (z). 

The St re-sums the threshold logarithms In 2 x appearing in the hard kernels to all orders and 
it has been parameterized as 

2 1+2 T(3/2 + c) 



S t (x) 



[x(i-x)r, 



(A14) 



v^Fr(l + c) 

with c = 0.4. In the nonfactorizable contributions, St(x) gives a very small numerical effect to the 
amplitude |47J] - Therefore, we drop S t (x) in h n and h na . 

The evolution factors E$ and Ea entering in the expressions for the matrix elements (see 
section 3) are given by 

E e (t) = a s (t)exp[-S B (t) - S 3 (t)], E' e (t) = a s (t) exp[-S B (t) - S 2 (t) - S 3 (t)]\ bl = b3 , (A15) 
E a (t) = a s (t) exp[-S 2 (t) - S 3 {t)\, E' a (t) = a s (t) exp[-S B (t) - S 2 (t) - S 3 (t)]\ b2=h , (A16) 



in which the Sudakov exponents are defined as 

5 f f dp, 



s B (t) 

S 2 (t) 



M Bs 



3 Ji/b L ^ 



s x 2 - 



s/2 



+ s[(l-x 2 )^-M 



l/6 2 V 



(A17) 
(A18) 



with the quark anomalous dimension 7 g = —a s /n. Replacing the kinematic variables of M 2 to M 3 
in 5-2, we can get the expression for S 3 . The explicit form for the function s(Q, b) is: 



s(Q,b) 



AW 



q In 



-q 



WA q - h ) + 



hx{2q) 



In (26) + 1 



+ 



Q _ 

h 

8/3? 



In 2 (2q) - In 2 (26) 



In 



hi 



(A19) 
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where the variables are defined by 



q = ln[Q/(V2A)}, S = ln[l/(6A)], 



(A20) 



and the coefficients A® and Pi are 




02 



153 - 19n 



24 



f 



(A21) 



Uf is the number of the quark flavors and 7^ is the Euler constant. We will use the one-loop 
running coupling constant, i.e. we pick up the four terms in the first line of the expression for the 
function s(Q,b). 

APPENDIX B: ANALYTIC FORMULAE FOR THE B s PP DECAY AMPLITUDES 

Before we give the analytic formulae for B s — > PP decays, we analyze the amplitudes in some 
special cases which can simplify the formulae. 

• In the non-factorizable emission diagrams, the formulae can be simplified if the emission 
meson M 2 is it, r\ or rj' '. Their distribution amplitudes <p A (x) and 4> p (x) are symmetric and 
4> T (x) is antisymmetric under the change x — > 1 — x. From Eq. (|54|) we can see that the 
amplitude is identically zero for the (V — A)(V + A) operators; From Eq. (153]) and Eq. 
f[5"5"j) we can see that the (V — A)(V — A) contribution is the same as the (S — P)(S + P) 
contribution which can be simplified as 



• In the annihilation factorizable diagrams, the (V — A)(V — A) and (V — A)(V + A) operators 
give the same contribution, both from the vector current. If the final state mesons are 
charge conjugate with each other, the (V — A)(V — A) and (V — A)(V + A) operators 
give identically zero contributions due to the conservation of the vector current. Such 
operators do not contribute to the B s — > Ktt decay either, in the 577(3) limit. The 577(3) 
symmetry breaking, i.e. the difference in the distribution amplitudes of tt and K meson, can 



M 



rLL 




x h n (xi, x 2 , £3, h, b 2 ) - x 3 03 (x 3 ) + 2r 3 x 3 4>J(x 3 ) ai(t' b )E' e (t' b ) 




(Bl) 
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induce small deviations from zero. We expect this kind of operators can not give important 
contributions to the branching ratios. 



• In the nonfactorizable annihilation amplitudes of B s — > nn, B s — > r)(nn)r)(nn) and 
B s — > r)(ss)r)(ss), the (V — A)(V — A) operators and (S — P)(S + P) operators give equal 
contributions as can be seen from the factorization formulae M ann by making X2 *-> 1 — £3. 

1. The case without rf 1 ' 



Tree operator dominant decays: 

G p 



V2 

-^V tb V; d { f n F^ K [04 + a 10 ] + UF§^ K [a 6 + a 8 ] 



+M^ K [C 3 + C g ]+f Bs F ( 



LL 
ann 



04 — 2 010 



+M, 



LL 



C^ — — Cq 

6 2 



+ M, 



LR 

ann 



Ck — — CV 
2 



fn i 71 



SP 

ann 



a 6 — 



(B2) 



71 



3 3 1 



-a 6 + -a 8 



ann 



-a 6 + -a 8 



+ fnF ( 



LL 
s*- ann 



3 1 

A 2 2 



-a 4 + -a 10 



LL 



-c 3 + -c 9 



M. 



LR 



-c 5 + -c 7 



-c 



10 



(B3) 



Pure annihilation type decays: 



v/2 



5l 

v/2 



2C4 -|- 2Cg 



:C 8 + -C- 



10 



■ (B4) 



51 



QCD penguin operator dominant decays: 



A(B° S -> K~K~ 



%K 6 TC{/^m [oil + M^ K [Ci] + M££[C a ] 

"71 



V^{/*Fg^ [04 + a 10 ] + f K F§?_+ K [a 6 + a 8 ] 



+M^ X [C 3 + C 9 ] + M|« [C 5 + C 7 ] + / fl .F t 



ann 



° 6 ~~ 2° 8 



LL 



M 



LP 



C 5 - -C 7 



+ M, 



SP 

ann 



+ {ML L n [Ca + c 10 ] + mZ [C 6 + c 8 ]) K _^ + }, 



(B5) 



+M 



Gf , 

LL 



04 — 2° 10 



f P 5P 



— —Co 



M LR 



Ck — — C7 



SP 
ann 



a 6 - 



+M, 



LL 



-Cg + C4 — 2^ 10 



M LR 

ann 



c 5 - ~c 7 



(B6) 



+M, 



LL 



1 



C4 — 2 < ^ 10 



+ M, 



rSP 



2. The case with rf'' 

As discussed in the last section, we use the quark flavor basis for the mixing of r\ and r\i. So 
we divide the amplitudes into the r\ n = (uu + dd)/\^2 and ss component. 



V2A(B° S ^ Vn K°) = ^V ub V: d {f n F^ K [a 2 ]+M^ K [C 2 ]} 



G 

72 
Gf 

71 



-- v tbV t * d ^f n F^_ >K 



a& - 2 a » 



fn F 



LL 

ann 



04 — 2" ai ° 



^LL 



+fnF B . K 



1 1 1 

2a 3 + a 4 - 2a 5 - -a 7 + -a 9 - -a 10 

1 1 1 

C3 + 2C4 + 2C*6 + -Cg — -Cg + 2 10 



+fs 3 F ( 



SP 

ann 



a& - 2°8 



LL 

ann 



— — Co 
2 



1 x nnn 



5 2 



•,(B7) 
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LL 

B S -,K 



+ fi<FB^ Vs 



ct>4 — 2 ai ° 



+fKF§^ Vs 



a 6 — 



+ f B F SP 

1 J &s ann 



C'\ — —Co 

6 2 



a 6 — 



+ M L B L 



C4 + Ce — ~C$ — 2^ 10 



5 2 



LL 

ann 



-aio 



LL 



C'? — —Co 
d 2 



M, 



LB 



Cr — — Cn 

5 2 



(B8) 



The decay amplitudes for the physical states are then 



A{B S -> r/A'°) = -> Vn K°) cos0 - 

A(B S -> v'K ) = A(B° S Vn K°) sin</> + A(B° S 



K°7] s ) sin 0, 

K°7] s ) cos (J). 



(B9) 
(BIO) 



For B s — > rj^TT the decay amplitudes are defined similarly, 



A(B S -> A) = -> An) cos, 

A(B S -> A') = A(B° S -> An) sin 



7r 77 S J sin<; 

7T°?7 S ) COS i 



(Bll) 
(B12) 



where 



■^^fe^is'i fB s Knn 



3 3 
2° 9 + 2° 7 



V 



Af. 



LL 



2 Cl ° 



,(B13) 



and 



V2A(B° S 



A. 



G p 

71 



-k V ubV* s { UF^ Vs [a 2 ] + M^ Vs [C 2 ] 



— CLq — — CLj 

2 2 



+ M LL 



B s —>Vs 



(B14) 
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For B s — > we have 



V2j4(£ s ° 77^, 



^Jc^ Us ann L 

^|K 6 K* S {/^IX [oa] + M| s % s [C 2 ] } 



LL 

B s -+r) s 



1 1 

2a 3 - 2a 5 - -a 7 + -a 9 



+M 



LL 



2C4 + 2C*6 + -C*8 + 2^ 10 



V2G F V tb V t :\f s F^ v 

+f s Fr 1 

-f R F s 



a 3 + «4 - ci5 + - ~ 2 ai ° 



SP 



a 6 — - Gs 



(B15) 



(B16) 



C3 + C4 + C*6 — — Cg — — C9 — — c 



10 



,S'L 



° 6 ~~ 2 a§ 



+ M, 



LL 

ann 



C3 + C4 + Cq — -Cg — -Cg — 2^1° 



117) 



The decay amplitudes for the physical states are then 



y/2A{B s - 


-> 7777) 


= A(B° S - 


A(B S - 


■+ 7777') 


= [A(B° S 


V2A(B S - 


► 7777) 


= A(B° S - 



?7 n ?7 n ) cos 2 + A(B° S -> 77,77,) sin 2 



* VnVn 

rj n rj n ) sin 



A(5° - 77,77, 
2 - + A(B° S - 



cos </> sm 1 
77,77,) cos 2 ( 



sin(20)A(B s ° 77 n 77,)(B18) 
77 n 77,) cos(20fB19) 
sm{2<j>)A{B° s -> 77 n 77,)(B20) 



APPENDIX C: ANALYTIC FORMULAE FOR THE 5, DECAY AMPLITUDES 



1. The tree dominant decays 



7T 



v/2 



— -^W^ UFbs-^k* [°4 + «io] - UFi^ K * [a 6 + a 8 ] 



+M|f_^ [C3 + C 9 ] + / Ss F c 



LL 

ann 



04 — 2" a io 



SP 

ann 



a 6 ~ 2°8 



+M, 



LL 



— — Cq 

d 2 



Lii 



C 5 - 2^ 



(CI) 
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A(B° S -> p-K- 



G p 
G p 



V tb V; d [f p F^ K [a 4 + a 10 ] + M#_ K [C 3 + C 9 



v/2 

+M|* [C 5 + C 7 ] + f Bs F ( 



LL 

ann 



a 4 — 2 ai o 



SP 
s-"- ann 



+M, 



LL 

ann 



Cq — — Cq 

6 2 



M. 



LR 



Cfi — C-! 

5 2 



(C2) 



V2A(5° -> 7r°A'*°) 



3 1 3 
— a 4 — -07 + 2" ai ° + 2 fl9 



+f B F L 
+M, 



LL 
s-"- ann 



-a 6 + -a 8 



-04 + 2" ai ° 



- fs F SP 

J Us ann 



3 13 
-C3 + + 2^9 + 2^ 10 



-06 + r ( 



LL 

ann 



"C3 + 2^9 



— M 



LL? 



-c 5 + -c 7 



(C3) 



C 



^ pjF 



3 1 3 
-a 4 + -a 7 + -aio + -a 9 



LL 

Bs^K 



+M 

1 J iJs^ ann 



1 3 

— —Co -f- — Cin 

d 2 J 2 



" a 6 + 2° 8 



M 



LL 



m B,^K 



"C3 + -Co 



:C 8 



LL 
s"- ann 



ff> F 



-C 6 + -C 7 



-a 4 + -aio 



+ M L 



ann 



-c 5 + -c 7 



(C4) 
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V2A(B° S ^ uK°) = ^V ub vdf w F^ K [a 2 ]+M^ K [C 2 ]} 



V2 



1 1 1 

2a 3 + a 4 + 2a 5 + -a 7 + -a 9 - -a 10 



+M 



LL 

Bs^K 



C*3 + 2C*4 — -Cg + 2^10 



LR 

B S ^K 



5 2 



-M 



2C*6 + 2*^8 



LL 

ann 



«4 — 2" ai ° 



ann 



a 6 — 2° 8 



+M, 



LL 



C3 — 2^ 9 



+ M, 



LP 



C*5 — — C7 



(C5) 



2. The pure annihilation type decays 



^ P ) = 



^V ub V: s U Bs F^ n [a 2 ] + M a L n L n [C 2 ]| - ^ Vtb V; s ^f Bs F a L n L n [a 3 + a 9 ] 
-/b,^ [a B + or] + M a L n L n [C 4 + C w ] - M s J n [C 6 + C 8 ] + [tt+ «-> p 



(C6) 



A(B° S -> p + 7r-) 



^K C| /b Ann [a 2 ] + M^[Cy 1 - / Bs C [a 8 + a 9 ] 



2A(£° -> 7T p°) 



-fB.F% [a 5 + or] + [C 4 + C 10 ] - [C 6 + C 8 ] + [p + <-> n 



A(B° S 
Gp . 



7T+p-) + A{B° S -> p + 7T 



10 



M. 



SP 



(C7) 
(C8) 
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3. The QCD penguin dominant decays 



-^V tb V*^f K F^ K * [a, + o 10 ] - f K Fl^ K , [oe + og] 
+M| s L _^ [C s + C 9 ] + / B F a L n L n 

+ fBsFann [ a 3 ~ «5 ~ «7 + O^] 



1 1 1 

a.3 + 04 — a 5 + 2 fl 7 ~ 2° 9 — 2 ai0 



-Af 



Li? 



1 



C 5 - -c 7 



-Mb%k* [C 5 + C 7 ] - f Bs F ( 
+M LL 

1 nnr, 



SP 
ann 



° 6 ~~ 2° 8 



-M 



Cr — — c« 



(CIO) 



V ub V: s {f K *F^ K [ ai ] + M^iCr] + / fl .i^ [a 2 ] + M^ n [C 2 



Gf 

-%V tb VZ{ faFgUx [a* + o 10 ] + [C s + C 9 ] + Mj£+ K [C B + C 

+ fpF LL 

' J ts ann 



03 + 04 — 05 + —07 — -Og — 2 Ol 



J £>s ann 



+fB.F™ [03 - a 5 - a 7 + anW* + M, 



Lii 



c 5 - -c 7 



a 6 — 2° 8 



r LL 
'ann 



C3 — -Cg + C4 — -jC\Q 



M, 



SP 



C*6 — ~C$ 



+M, 

[^4 + C 10 } K ,^ K - M^ n [C 6 + C 8 j^,^^ ?, 



(Cll) 
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71 



VMf K F™ 



04 — -oio 



a 6 — -08 



+M 



LL 



-Mil*, 



c 3 - 




+ f R A iL 

1 J Ms^ ann 




K 


- fs F SP 

J £>s ann 



1 1 1 

a 3 + a 4 - a 5 + -a 7 - ~a 9 - -a w 



o& ~ 



LL 

ann 



M 



LL 



C 4 — -C ho 



-M. 



LR 



c 5 - \c 7 



M 



SP 



C*6 — 



[if* <-> A'] 



LL 

ann 



1 1 

^3 - ^5 + 2 fl 7 - 2° 9 



(C12) 



G 

'71 



04 — -aio 



C^ — — Co 



+M LR 



Bs^K 



+f Bs F c 
+f B F c 



SP 
s " ann 



c 5 - -c 7 



06 -gt 



+ f R A LL 

1 J iJs^ ann 

+ M, 



a 3 + a 4 - a 5 + -a 7 - -a 9 - -ai 



fLL 
"ann 



1 1 

C*3 — -Cg + C 4 — 2^10 



Lii 

ann 



c 5 - -c 7 



M 



SP 



C% — -Cs 



[K* <-> A'] 



LL 

s ann 



0-3-0-5+ 2°7 - 2 fl9 



LL 



K*^K 



C*4 — 2^ 10 



(C13) 



1 1 

«3 + a 5 - ~a 7 - -ag 



LL 



SP 



^6 - 2° 8 



M 



LL 

Bs^K 



Ca — l^Cio 



+ fx-Fft 



04, — 2" ai o 



C, 



:Co 



-M SP 

Bs^fK 



-fn A 



SP 

ann 



a 6 



2 
1 

— i 
2 



:C, 



-c 7 



+ M, 



LL 



C*3 — -Cg 



fn A 

M L 

x VJ - n.i 



LL 
s *■ ann 



Oi — 2" ai ° 



tLR 
*-ann 



C 5 - -Cr 



(C14) 
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4. The Electroweak penguin dominant decays 



G 



Li 



-Og GS7 



Mjf_ 



:C 8 + -C 



10 



5. The decays involving 7] and 77' 



V^A^^^iT ) = ^V ub V: d {f n F^ K ,[a 2 } + M L B ^ K ,[C 2 ]) 



V2 
Gp 



1 1 1 

2a 3 + a 4 - 2a 5 - -a 7 + -a 9 - -ai 



C3 + 2C*4 — — Cg + 2 < ^ 10 



2Ck + — Co 



LL 

nn 



1 

^4 — 2° 10 



f R F 



SP 
3" ann 



+M, 



G 



LL 



C'? — — Cq 

d 2 y 



C 5 - 2^ 



a 3 -a 5 + ~a 7 - -a 9 



+ fK*Fg^ n 



a 1 



+M 



LL 

Bs^K* 



+M SP 



B S ^K* 



C4 — -C '10 



C*6 — -Cg 



-Ds >T}S 



C3 — —Cq 

6 2 



C5 — — C7 

a 2 



LL 

nn 



a 1 



+fB s F ( 



SP 
ann 



° 6 _ 2° 8 



LL 



C3 — ~Cg 



LR 



C 5 - -C 7 



The decay amplitudes for the physical states are then 

A{B S -> r/A'*°) = A(B S ° -> r/ n iT°) cos0 - A(B° S -> r/ s iT°) sin, 
-> 7]'K*°) = A{B° S -> r/ n iT°) sin</> + A(E S ° -> r/ s iT°) cos, 

For 5 S — > r/'^p the decay amplitudes are defined similarly, 

p\) = -> pV) cos0 - -> pV) sin0, 

A(B S -> p°rj') = A(B° S - pV) sin0 + -> p\ s ) cos0, 
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where 



2A(B° S - An) = f B F c 



LL 

ann 



2 ag " 2° 7 



M, 



LL 



2 Cl ° 



3 

-( 

2 



G, 



and 



V2A(B° S -> p%) 



G 



^fW^ M + Mjf_^ [G 2 ] 



3 3 
-«7 + 2 fl 9 



+ M 



LL 



2 Cl ° 



For decays 5 S — > rj^'u, we have 



Gf 
72 



/i^ann M + M^ n [G 2 ] 



LL 



2C4 + 2^!° 



— M 



SP 



+ -Gg 



+/p s ^ ( 



LL 
ann 



1 1 

2a 3 - 2a 5 - -a 7 + -a 9 



[?? n <-> w] 



G 



V^A(S? - = ^§K 6 C j [a 2 ] + M^ Vs [G 2 ] 



2a 3 + 2a 5 + ^a 7 + ^a 9 



+M 



LL 

B s —Hf]s 



2G4 + gGlO 



^1% 



2G 6 + -G 8 



And the decay amplitudes for physical states are 



A(B a -> 77^) = -> 77^) cos0 - A(B° S -> cu?7 s ) sin0, 

— > UT] S ) cos0. 
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For decays B s — > -qcj), we have 



C 



V2A(B° S -> r) n <f>) = ^V ub V:A f n F^ [a 2 ] + M| s % [C 2 ] 



G^ 



1 1 

2a 3 - 2a 5 - -a 7 + -a 9 



LL 

B a - 



2C4 + 2^10 



M 



SP 

B s - 



2C e + -C 8 



G f 



1 1 1 

(13 + 04 — a 5 + -07 — -ag — 2" ai o 



+M S / 



l B a 



a 6 — 



Cfi — —Co 



M; 



LL 

B s - 



+ fB s F c 



LL 

s arm 



+M, 



LL 



C3 + C 4 — -Cg — ^GlO 



SP 
ann 



-M 



LP 



c 5 - -c 7 



- M 



SP 



Gs — -Gs 



a 6 — 



fas 



(C28) 



C3 + C4 — -Cg — ^Gio 



1 1 1 

a 3 + a 4 - a 5 + -a 7 - -a 9 - -a w 



(C29) 



And the decay amplitudes for physical states are 



A(B S -> 770) = A(5 S ° -> r? n 0) cos0 - -> rj s (j)) sin0, 
A(5 S -> r/0) = A(B° S -> r/ n 0) sin0 + ,4(5° -> r/ s 0) cos</>. 



(C30) 
(C31) 



APPENDIX D: ANALYTIC FORMULAE FOR THE B s -► DECAY AMPLITUDES 

There are three kinds of polarizations in the B s meson decays to two vector final states, namely: 
Longitudinal (L), parallel (N) and transverse (T). The decay amplitudes are classified accordingly, 
with i = L,N, T. 
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1. Tree dominant decays 



-^VtbVtd^fpF^^ [a 4 + a w ] 
+M£i K * [C 3 + C 9 ] - M L B fi K * [C 5 + C r ] 



i J as*- ann 

+M, 



0,4 — 2°!° 



J £>s ann 



a 6 — 



rLL,i 
'ann 



d 2 



M 



5 2 



(Dl) 



f P F^ K , [a 2 ] + M^[C 2 ] 



C^ 



LLA 



V2 

+M L B L ' K , 



3 3 1 
-a 4 + -a 7 + -a 9 + -ai 



1 3 

-C3 + -Cg + 2^!0 



M 



3 c 

2^ 8 



J r>s c 



LL,i 
ann 



-CL4 + 2" ai o 



ann 



1 

-a 6 + 2 ( 



+M, 



-C 3 + -Cg 



-c 5 + -c 7 



M 



LR,i 



-C 5 + 2^ 



(D2) 



*0\ 



G 
72 

G p 

+M LJ 



V tbV t *A fuF B ^ K * 



1 1 1 

2a 3 + a 4 + 2a 5 + -a 7 + -a 9 - -a w 



-M 



B S ^K* 

LR,i 
B„->K* 



C3 + 2C*4 — —Cg + ^^iO 



- f R F LL ' { 

J £>s ann 



c 5 - -c 7 



a 4 — 2" ai ° 



M 



SPA 
B S ^K* 



fB.F t 



SP,i 
ann 



2C 6 + -Cg 



a 6 - 2 ( 



+Af 



C? — — Cq 

2 



M. 



LRA 



C 5 - -C 7 



(D3) 
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Pure annihilation type decays 



-^V tb V*if B F c 



LL,i 
s-"- ann 



2a 3 + -a 9 



fn F 

J 13 s c 



LR,i 
ann 



2a 5 + 



+M, 



LLA 



2C4 + -C \q 



+ M, 



SP,i 
ann 



2C 6 + -C 8 



V2A l (B' s 



P°P°) 



Gp 
71 



VubV*A f B FL L n l N + M^[C 2 ] 



LLA 



--^=VtbVt s { fB 3 F ann 



2a 3 + -a g 



fp> F 



LR,i 
ann 



2a 5 + 



LLA 



2C4 + 2^10 



2C*6 + 2^8 



H = ^K&K^ /aFS* [02] + M££[C 2 ] 



_ — l ^tfe^ts^ fBs F ann 



2a 3 + -a 9 



+M, 



LL,i 
ann 



2C4 + 2^10 



1 x nnn 



fsF ( 
2Cfi + — Cg 



LP,i 
ann 



2a 5 + -a 7 



A) 



r 1 Gp 

-^^ubV* s <yfB a Fann M + M^ l [C 2 ]j ~ -j=V tb V t * s \ fsF^ 



rag 



+fs a F ( 



LR,i 
s" ann 



-07 



ann 



-c 



10 



M. 



SP,i 



:C 8 
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QCD penguin dominant decays 



A*(B° -> K*~K* + ) 



^V ub V; is {f K *F L B ^ K , [ ai ] + Mg^jf, [Ci] + fB 3 Fann N 



+M^/[C 2 ] f - /j^i^i. [04 + a 10 ] 



[C 3 + Cfe] - [C B + C 



LR,i 



+ f r F LL,i 

1 J £>s A arm 
■/ t>s ann 



1 1 

a 3 + a4 — -a 9 — 2 a io 



a 6 



+ M, 



LL,i 
ann 



fs s F ( 
1 



LR,t 
s~ ann 



«5 - 



Co + Cm 



-C 



10 



_ M LR,i 



c 5 - -c 7 



M, 



SPA 



c fi 



+ (/i^S* [flS + Of,] + [«5 + « 7 ] 



[C 4 + C; 



10 



5P,i 



[C6 + C 8 ])^ 
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?LL,i 
B.^K* 



1 

04 — 2 Ql ° 



-r ivi B Rt 



C 3 — — Cg 

d 2 



m B s ->K* 



-fn F 

J 13 s £ 



ann 



c 5 - -c 7 



M LL,i 



LL.i 

s ann 



1 1 

a 3 + a 4 - -a 9 - -a 10 



1 1 



-M LR,i 

ann 



C5 — — C7 

5 2 



ann 



1 Jj>! ann 

+ IM 



fB.F t 



LR.i 

ann 



a-b - 2°7 



fLL,i 
"ann 



Ca — —C 10 


1 ann 


c 6 


-K 















(D9) 



V2' 



1 1 

«3 + «5 - 2°7 - 2°9 



-l-M LL,i 



Ca — 2^ 10 



Ck — ~C-[ 

5 2 



r LL,i 



J £>3 ann 



C? — — Cq 
d 2 9 



04 — 2" ai ° 



p rpLL,i 



B a —*K* 



a a — 2 010 



C« — —Co 



- f B F SP ' { 

J £>s ann 



a& - -a 8 



C 3 — — Cg 
d 2 y 



M. 



C 5 - -C 7 



(D10) 



C^ 
71 



/u^iiL [02] + Mj£L [C 2 ] 



LLy 



-^v tb v; s \ UF^ 



1 1 

2a 3 + 2a 5 + -a 7 + -a 9 



2C 4 + 



M 



SPa 



2C6 + -Cs 



(Dll) 



2Gf 

'7f 



LL,i 



1 1 1 

a 3 + a 4 + a 5 - -a 7 - -a 9 - -ai 



+M 



-M 



LL,i 
B s ^ 

SP,i 
B a ^<, 



C 3 + Ca — -Cg — 2 < ^ 10 



M 



LR,i 
B s ^d 



Ck — — C7 

& 2 



C6 — 



x pLR,i 
J -Os ann 



as - 2 a? 



+ fn F ' 

1 J £>s ann 

- fs F SF 

J Us am 



1 1 

a 3 + 04 — -a 9 — 2" ai o 



1 

a 6 — 2 a§ 



— M 



c 5 - -c 7 



1 1 

C 3 + C 4 — -Cg — ^CiQ 



Cfi — —Co 



(D12) 
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4. Electroweak penguin dominant decays 



y/2A\B° a -> p\ 



Gp 

71 



Vu b V:Af p F^M+M^,[C 2 



LLA 



(a 9 + a 7 ) 



I/I/,' 



-c 



10 



M 



(D13) 
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